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ABSTRACT 


A long standing concern of the Navy has been the need to reduce 
structural vibrations in plates, such as ship's hulls. Recently, it was 
proposed to use waveguide absorbers as a means to reduce structural 
vibrations. A waveguide absorber is a device which is made from damping 
material but mounted to the structure at a point rather than the full 
surface. A waveguide absorber removes and then dissipates vibrational 
energy from the structure through traveling waves. The performance of a 
waveguide absorber is determined from the driving point impedance, the 
ratio of force to velocity at the attachment point. This study has determined 
the theoretical driving point impedance for viscoelastic circular plates. 

The constitutive relations for viscoelastic material were evaluated 
and it was proven that a complex Young's modulus and shear modulus 
can be used to represent the viscoelastic material. Mindlin's theory for 
elastic isotropic plates, with complex moduli, was used to solve for the 
driving point impedance. The solutions are in terms of Bessel functions. 
Poisson's boundary conditions were used at the free edge of the circular 
plate, and clamped boundary conditions were used at the attachment point, 
which is at the center of the plate. A non-numeric language, REDUCE, 
was used to solve for the driving point impedance given the appropriate 
equations of motion and boundary conditions. The complex Bessel 
functions, needed for the results of the REDUCE program, were generated 
by a FORTRAN program. To validate the results of this study, elastic and 
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viscoelastic plates were tested to determine their driving point impedance. 
A comparison of the theoretical results and the experimental results shows 


that there is agreement for all cases studied. 
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I. INTRODUCTION 


A. BACKGROUND 

The reduction of structural vibration in ships and submarines has 
been a long standing concern of the Navy. Reducing structural vibrations 
can increase the life of equipment as well as decrease the amount radiated 
noise into the water. One of the most widely used methods of passive 
structural vibration reduction has been the addition of damping material 
directly to the surface of the structure[Ref. 1 : pp. 312 -362]. Recently Ungar 
and Kurzweil [ Ref. 2 ] have proposed the use of a waveguide absorber to 
reduce structural vibrations. A waveguide absorber is a damping device 
that is attached indirectly to the structure through a mount, at a point, 
rather than directly to the surface. The purpose of the waveguide absorber 
is to remove energy from the structure and then dissipate the energy as 
traveling waves. The waveguide absorber is expected to remove more 
energy per unit mass of damping material than a surface damping 
treatment can. 

Past studies of waveguide absorbers have been by Ungar and 
Kurzweil [Ref. 3] , by Ungar and Williams [Ref.4], and by Lee [Ref. 5] . 
These studies show that a waveguide absorber is able to reduce the 
magnitude of the vibrations in a structure. The theoretical results are 
based upon the driving point impedance, the ratio of force to velocity at the 
attachment point. There has not been good agreement between the 


experimental results and the theoretical results in these studies. However, 


in the study by Lee [Ref. 5] of viscoelastic beams, it was determined that 
using an underlining theory that included the shear deformation and 
rotatory inertia, i.e. the Timoshenko beam theory, produced better results 


than those theories which did not include them, i.e., the Euler beam theory. 


B. PURPOSE 

Based upon the experimental success of beam waveguide absorbers to 
reduce structural vibrations, there is reason to believe that a two 
dimensional geometry may provide increased damping per unit mass of 
damping material than the one dimensional beam does. The purpose of 
this study is to investigate the feasibility of using viscoelastic circular plates 
as waveguide absorbers, by determining the driving point impedance of the 
plate. The driving point impedance will be evaluated over a frequency 
range of 10 to 2000 Hz. This range of frequencies was chosen because data 
gathering techniques and material information are not reliable below five to 
ten hertz, and the distance that vibrations or noise can travel in water above 
2000 Hz is very small and normally ignored. It is the goal of this study that 
the driving point impedance be presented an easy to use equation, or an 
easy to use computer program. The driving point impedance should be able 


to be calculated with a limited number of inputs. 


C. FORMAT 
This study has been broken down into four sections even though the 
actions occurred at the same time. The sections are theoretical results, 


experimental results, a comparison of the theoretical and experimental 


results, with conclusions, and recommendations for future studies. The 
body of this study was written based upon all of the plates used, but only 
shows, in the figures, the results for one plate size. The figures for the 
other plate sizes are found in Appendix C. This was done so as to keep the 


body of this study as small as possible. 


IL THEORY 


A. VISCOELASTIC MATERIAL MODELING 

In the past study by Bailey and Chen [Ref. 6], it was shown that one of 
the best methods of modeling viscoelastic material is by separating the 
elastic and viscoelastic parts. The elastic portion is not time dependent, 
where as the viscoelastic portion is time dependent. The model of 
viscoelastic material can be created by using a resistive element as shown 
in Figure 1 and a damping element as shown in Figure 2, where o is the 
stress, € is the strain, E is the Young's modulus, and C is the damping of 
the material. These models are based upon electrical impedance where 
strain is the equivalent of voltage, stress the equivalent of current, and the 


inverse of the characteristics is the equivalent of impedance. 
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Figure L Resistive Element Model. 
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Figure 2. Damping Element Model 


The stress versus strain relationship for the resistance elements is: 


=o 
c“f (2.1) 
and for the damping element is: 
_|o 
a ig (2.2) 


The combination of a resistive element and a seonies element is known as 
the Maxwell model [Ref. 7]. The model shown in Figure 3 is based on a two 
element Maxwell model, where the full element, with damping, represents 
the viscoelastic portion and the degenerate element, no damping, 
represents the elastic portion of the material. The notation used is, Ee is 
the elastic Young's modulus, Ey is the viscoelastic Young's modulus, Cy is 
the viscoelastic damping, 6e is the elastic stress, Oy is the viscoelastic 


stress, 01s the total stress, and € is the total strain in the material. 








Figure 3. Viscoelastic Material Model. 
The total stress,o, of the material is: 
Oe ee (2.3) 


The strain versus stress relationship for the elastic portion of the matenal 


is: 


6} 
a cele 


or in a more familiar form: 


O,=He (2.4b) 


For the viscoelastic portion the strain verse stress relationship is: 


mela le 
E= E. + in (2.5a) 


or taking the time derivative of equation (2.5a) gives: 


O Oo 
V V 


It is concluded that the total stress of equation (2.3) can be rewritten as: 


Oo=E,e€+ 0, (2.6) 
where the viscoelastic stress must obey the relationship of equation (2.5b). 
The stress versus strain relationship can be characterized as the normal 
elastic relationship plus an added term that incorporates the viscoelastic 
effects. 

For reference purposes the coordinate system of Figure 4 is used for 


the discussion that follows. 





Figure 4. Rectangular Coordinate System for a Plate. 


For thin plates, it is normal to assume that the plate undergoes a biaxial 
stress state [Ref. 8], but for this study triaxial stress is used with the 
modification that the strain in the z direction is zero and that the stress in 
the z direction is very small and can be ignored. This assumption will 
allow the inclusions of the shear deformation and rotatory inertia later. 


The stresses acting a small piece of material are shown in Figure 5, where 


Oj, Or 0}j, is the normal stress on the 1 face in the i direction, t t is the shear 


stress on the i face in the j direction, which is equal to t ji- 





Figure 5. Stresses in a Cube of Material. 


Based upon the model in Figure 3 with resultant equation (2.6), the three 
dimensional stress versus strain relationships for isotropic material 


become: 


E, 
o,= (ExtV ely) + Ox, 
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(2.7a) 
(2.7b) 


(2.7c) 


(2.7d) 


(2.76) 


(2.7f) 


where €j is the strain in the i direction, ¥ ij is the rotation in the i-j plane, 


which is equal to y jp and the sub-sub-script of v is for the viscoelastic 


portion. The viscoelastic terms must also meet the following conditions: 
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(2.8a) 


(2.8b) 


(2.8c) 


(2.8d) 


(2.8e) 


The use of these constitutive relations proves to be quite cumbersome 
in all but the most simplest of problems. This is because simultaneous 
differential equations must be solved to relate the stresses to strains. By the 
use of a complex model, where the time dependent terms are replaced by 
complex terms, some of these problems can be overcome. Figure 6 shows 


the Laplace Transform [Ref. 9], of the model shown in Figure 3. 
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Figure 6. Laplace Transform of Viscoelastic Material Model 


By applying circuit analysis methods to the model shown in Figure 6, an 
equivalent model, shown in Figure 7, with a complex Young's modulus, 
E*, can be determined. The viscoelastic elements are in series with each 
other and together are in parallel with the elastic element. Solving for E* 


yields: 
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1 1 
=a (f+4) (2.9) 


m E, \E, sCy 


aus 
1 








e 
Rearranging terms and substituting jw for s, where j= V- 1 and wis the 


rotational frequency in rads/sec, in equation (2.9) yields: 


2 24 22 2 

» EE, +o CLE +@ Cae oC EY (2.10) 
~ 2 92 V3 9,2 
as Oe Pe ones 


The resultant complex model is shown in Figure 7. 





Figure 7. Complex Viscoelastic Material Model. 


The complex model not only allows the use all of the elastic theory, with its 
well defined mathematics, but has the added benefit that most of the 
manufactures of viscoelastic material provide material data in this form 
(Ref. 1, pp. 363-440]. Instead of leaving the complex Young's modulus in 
terms of the elastic and viscoelastic parameters, as in equation (2.10), three 


new terms are defined and used. The real part, E’, is defined as: 


2 Dias 9 2 
= He + @ CE. + CLE, (2.11) 


Z 9 2 
Pe once 


ae 


The imaginary part, E", by: 
Ke” oC EK. (27) 


The ratio of the imaginary part to the real part is known as the loss 
damping factor, n, and defined as: 


oe 
= (2.153) 
E 
The complex Young's modulus may now be written as: 
E =E+jE (2.14) 
or as: 
E = E(1 + jn) (2.15) 


It is normal to assume that the Poisson's ratio,v , for the viscoelastic 
material is the same for the elastic and viscoelastic portions, which leads to 
the result that the shear modulus is also a complex parameter, G*, and is 


related to the complex Young's modulus by: 


x 


* E 


G 


The real and imaginary parts, respectively, of the Young's modulus and 
shear modulus are related in identical manner. The loss factor is the same 


for both moduli. 
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Elastic theory can now be used to analyze structures using 
viscoelastic material when Young's modulus and shear modulus are 
replaced by the viscoelastic's complex moduli. This is a great benefit in the 
analysis of structures with viscoelastic plates because a separate and new 


governing theory does not have to be created. 


B. ELASTIC SOLUTION 

For this study an elastic thin plate theory that includes the effects of 
shear deformation and rotatory inertia was needed to determine the driving 
point impedance of the viscoelastic plate. The use of such a theory, for 
beams, in the study by Lee [Ref. 5] provided better agreement between 
theoretical and experimental results than did theories which did not 
incorporate these effects. Mindlin [Ref. 10] has presented such a theory for 
isotropic thin plates. This work showed that isotropic elastic thin plates 
can be described by a set of three simultaneous partial differential 
equations. It is based upon assuming that the displacements in the x 
direction, u, and in the y direction, v, of Figure 4, vary linearly in the z 
direction and the displacement in the z direction, w, is not a function of z, 


or written in equation form: 


u=z-: W(x, y, t) (2.17a) 
v=z-Wy(x, y, t) (2.17b) 
w = W(x, y, t) (2.17c) 
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where wx and wy are potential functions to be determined. This implies 
that the displacements in the x and y directions, u and v, vary as shown in 


Figure 8. 





Plate Displacement 


Figure 8. Variations in x and y displacements as a function of z 
for the Mindlin Plate Theory. 


Mindlin's original work was in a rectangular coordinate system, but 
was adapted to cylindrical coordinate system by Mindlin and 
Deresiewicz [Ref. 11] and then updated by Irie, Yamada, and 
Aomura [Ref.12], and to curvilinear by Callahan [Ref.13]. The notation of 
Irie, Yamada, and Aomura [Ref. 12] is used in this study owing to studying 
circular plates and their use of dimensionless parameters. The reference 


directions of the cylindrical coordinate system are shown in Figure 9. 
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Figure 9. Cylindrical Coordinate System for a Plate. 


The displacements, in terms of cylindrical coordinates, which are 


equivalent to those of equation (2.17) are the radial displacement, ur, the 


angular displacement, ug, and the displacement in the z direction, w. 


They are expressed mathematically as: 


u,=z- Ww (r, 8, t) (2.18a) 
U, = 7Z- nes 6, t) (2.18b) 
w=W (r, @, t) (2.18c) 


where wrt and Wg; are potential functions that are to be determined. The 


moments and shears in a plate are shown in Figure 10 and defined in 


equation (2.19) [Ref. 10] for the case of cylindrical coordinate system. 
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M 9 
Figure 10. Moments and Shears in a Plate for Cylindrical 


Coordinates. 


(2.19a) 
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where Mjj, or Mj, is the bending moment on the i face per unit of length in 
the r-@ plane, Mjj, which is equal to Mjj, is the twisting moment on the 1 
face per unit of length in the r-@ plane, and Qj, is the shear on the i face per 


unit of length in the r-@ plane . The plate geometry is shown in Figure 8 


and Figure 9. 





Figure ll. Plan and Side View of Plate Showing Geometry. 


For the case of an isotropic plate in cylindrical coordinates, the 


Mindlin governing equations are [Ref.11]: 
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go ae 52 (2.20c) 


where p is the density of the material. It is not desirable to solve 
equations (2.20) in the present form, because of the simultaneous 
differential equations, but rather in a simpler form. Mindlin [Ref. 10] 
showed that the three simultaneous second order partial differential 
equations, can be combined to produce a single fourth order partial 
differential equation, in terms of w. The only drawback of using this single 
equation is that it can not handle the full number of geometric boundary 
conditions. Another suggestion by Mindlin [Ref. 10] was to assume that the 
displacements were all periodic functions of time of the type exp(ot). By 
doing this, the three simultaneous equations can be converted into three 
independent second order partial differential equations, which are much 
easier to solve, but still use of all the geometric boundary conditions. 

For this study, it was desired to use the geometric boundary 


conditions directly, in order to improve the agreement between the 
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theoretical and experimental results, and since most shipboard structural 
vibrations are from rotating equipment, periodic in nature, the final 
suggestion by Mindlin will be used. As such, the displacements are 


assumed to be periodic of the form exp(mt ), and expressed as: 


aot 
u_=z-W_(r, ee” (2.21a) 
jot 
ine oc V g(t O)e (221i) 
w = W(r, ge) (2.21¢) 


where the unknowns are Wr, We, and W. It is convenient to define the 


following dimensionless parameters that are based upon the plate geometry 


and the frequency [Ref. 12]: 
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R= (2) (2.22a) 
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RA - S 
where D is the flexural rigidity expressed by: 
3 
De =e (2.23) 
12(1 - v2) 


and x is the shear coefficient taken from Mindlin [Ref.10]. The shear 
coefficient is used to take into account that the shear is not a constant value 
as assumed, but varies across the thickness of the plate. For elastic plates 


the shear coefficient is normally set equal to: 


2 


= 2) (2.24) 


For the presently assumed case of periodic motion of isotropic plates, 


the governing partial differential equations become [Ref. 12]: 


2 2 

(Vv + 5. \w, = 0 (2.25a) 
2 2 

(v + 5. \w, =0 (2.25b) 
2 2 

(v +6 a =0 (2.25c) 


where V2 is the Laplacian operator and the variable w}, w2, and w3 are 


dummy variables which are related to the unknowns, yr, Wg, and W by the 


following relationships [Ref. 12]: 





Ow ow, 1 OW . see 
os (Ci ae ae ~+(6, ip ie Tr 90 ( : a) 
Ow Ow Ow 
_ eon Re 3 2.26b 
We=(% 1) 06 +(5,-1)z 06 .~=o or 
Was Bae (2.26c) 


It is much easier to work with the governing equations as shown in 
Equation (2.25) than those of Equation (2.20). The problem has been reduced 
from having to solve simultaneous partial differential equations, to solving 
three independent partial differential equations. Equation (2.25) is a form of 
Bessel's partial differential equation which has a solution in terms of 
Bessel functions [Ref. 14] [Ref.15, pp. 22-26]. For this study the solution to 


equation (2.25) is assumed to be: 


w,=[A,IJn( 8,4) +4 ¥n( 8,4) fos (n9) (2.27a) 
wo={AjJn( 8,2) +4 5¥n( 854) feos (n8) (2.27b) 
we={AJ (8.5) +4 .¥n( 85a) pin (n) (2.27c) 


where Aj are constants to be set by the boundary conditions, Jn(x) is the 
Bessel function of the first kind of order n, and Ypn(x) is the Bessel function 
of the second kind of order n. The Bessel function of the second kind of 
order n, has been added to the solution of Irie, Yamada, and Aomura 
[Ref.12] because of the hole in the center of the plates being used in this 
study. 

Six boundary conditions are needed to completely establish the 
solution to this plate problem. For the plates under study, three boundary 
conditions are at the center of the plate, at r=b, and three boundary 
conditions are at the edge of the plate, at r=a. The boundary conditions are 


given in terms of displacements and the stress in the plate. At the center, 
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r=b, for a plate being used as a waveguide absorber, it must be mounted to a 
fixture. It has been assumed that the fixture is rigid, has a radius of b, and 
clamps the plate such that there is no twisting and that the displacement in 
the z direction at r=b is the same for the fixture as the plate. In 


mathematical terms, the boundary conditions at the plate's center are: 


y ,(b, 8) =0 (2.28a) 
Y (b, 8) =0 (2.28b) 
W(b, 8) = Wo (2.28c) 


If the plate was not attached to the fixture, or if b=0, then the constants Ag, 
As, and Ag would have to be set to zero to insure that the solution remained 
bounded, as the Bessel functions of the second kind are unbounded as the 
argument approaches zero. 

At the free edge of the plate, r=a, the boundary conditions are no 
bending moments, twisting moments, or shearing on the free edge, or in 


mathematical terms: 


M _,(a, 8) =0 (2.29a) 
M ,,(a, 8) =0 (2.29b) 
Q (a, 0) = 0 (2.29c) 


These are known as Poisson's boundary conditions [Ref.16]. Poisson's 
boundary conditions can be used because the governing differential 
equations were not reduced to form a single differential equation which 
ignores, or reduces the effects of the shear deformation and rotary inertia. 
Had a reduced differential equation been used, then Kirchoff's 

boundary [Ref. 16] [Ref.17] conditions would have to be used. Kirchoff's 


boundary conditions for the free edge are: 


M ,(a, 8) =0 (2.30a) 


aM (r, 0 
eRe ae.) 4 _ (2.30b) 


Cre 8) +r 96 


r=a 


Kirchoffs boundary conditions are an attempt to include the three 
geometric boundary conditions of Poisson that exists at the free edge into 
two equations. The necessity for this arose from a plate theory, equivalent 
to the Euler beam theory, which could not incorporate three edge boundary 
conditions but only two. The present theory, Mindlin's, does not give up 
boundary conditions for simplicity. The use of a reduced theory and 
Kirchoffs boundary condition can lead to a large disparity between the 
theoretical and the experimental displacements [Ref. 18] when r 


approaches the free edge at a. 


All that is left to allow analyzing the plate is to relate the 
displacements, or potential functions, to the moments and shears. The 


equations below accomplish this [Ref. 12] : 








M,. = pie t vie a “o (2.31a) 
M oo = D1 Hy, + et) madd (2.31b) 

D 1( Vs Ne 
Mg IY, {2 36 -v,)* i — 
Q,= «Gh y, i: oN ) (2.31d) 
Qy= ch(y, ; io ) (2.31e) 


Finding the constants Aj to AG, is a straight forward matter of 
combining equation (2.27) into equation (2.26), and then into equation (2.31) 
from which the boundary conditions of equations (2.28) and (2.29) can be 
satisfied. Once the constants Aj to Ag are known, the potential 


functions, Wr and Yo, the displacements,uy, ug, and w, and the moments 


and shears, Mrr, Mogg, Mro, Qr, and Qp, are known anywhere in the plate. 


With this information the driving point impedance can be calculated. 


Mechanical impedance, Z, is defined as: 


F 
7 (2.32) 


where F is the force at the point of interest and V is the velocity at the point 
of interest. For this study, the point of interest is the the driving point 
where the fixture mounts to the plate. Based upon the assumption made 
about the fixture and nature of the displacements, the force, F, at the 
driving point can be calculated by: 


2n 


jot 
F= | 1Q,(b,0)e” 20 (2.33) 
0 


and the velocity, V, is related to the displacement at the driving point by: 


jot 
-. AV Gs )_ 


. ae 
ja W (b, ge?” =jow .e (2.34) 
Combining equations (2.33) and (2.34) into equation (2.32) yields: 


2r 


[rQ .(b, 6) 08 
_ 0 (2.38) 


Equation (2.35) is the driving point impedance. 
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The steps necessary to solve for the driving point impedance have 
now been laid out. The only difficulty lies in determining the constants Aj 
to Ag for any plate geometry (a,b,h) and material (p,E”,v). The work 
necessary to determine these constants is beyond the capability of normal 
hand written calculations. Assuming that the driving point impedance 
can be derived, the driving point impedance must be calculated over and 
over for different frequencies, material, and geometry. The only way to 
solve for the driving point impedance is with the aid of a computer. The use 
of the computer was split into two parts. The first was to use a symbolic 
language to algebraically determine the driving point impedance, and the 
second was to use a numerical language to calculate the driving point 


impedance for a given set of conditions. 


C. COMPUTER IMPLEMENTATION 

Figure 12 shows the flow path used in the symbolic program to 
determine the driving point impedance. One of the primary considerations 
in the selection of a non numerical language, which will be used to 
determine the driving point impedance, was the ease of differentiating 
Bessel functions. The language of REDUCE [Ref. 19] has the capability to be 


easily taught how differentiate functions, so was used in this study. 
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Methods for Differentiation 
of Bessel Functions 


Assumed Solution 
Eq (2.27) 


Determine Potential 
Functions and w 
Eq (2.26) 


Determine Moments 
and Shears 
Eq (2.31) 


Determine Boundary 
Condition Equations 
Eq (2.28) and (2.29) 


Solve for Al to A6 


Determine Driving 
Point Impedance, Z 
Eq (2.35) 


Write Z in FORTRAN 
Type Statements 


Figure 12. Symbolic Language Program Flow Chart. 


The REDUCE program did the differentiation of the Bessel functions 
and substitutions as required very well. The program showed that the 
resultant motion was symmetric with n equal to zero. This agrees with the 


boundary conditions imposed at the center of the plate, motion only in the z 


direction with no twisting. It was only at the point of trying to determine 
the constants Aj to Ag did REDUCE have trouble. In effect, the REDUCE 
program was solving a linear six by six symbolic matrix. In the present 
implementation of REDUCE on the Naval Postgraduate School's IBM 3033 
computer there was insufficient storage to to find all six constants directly. 
In order to alleviate this problem, the program was broken up into two 
parts, one part to set up the necessary equations to be solved, and the other 
to solve the simultaneous linear equations. Cramer's rule of linear matrix 
solution [Ref. 20] was employed to solve the simultaneous linear equations 
through the use of dummy variables. Even using this approach there was 
insufficient memory to solve the system of equations. The symmetry of the 
problem was then taken advantage of, A3 and Ag do not need to be solved 
because sin(n@) is zero. The result was, that the REDUCE program only 
had to solve a four by four matrix for Aj, A2, A4, and A5, which it was able 
to do. Once these constants were found then the magnitude of the force, F, 
magnitude of the velocity, V, and the driving point impedance, Z, were 
found. The driving point impedance is in terms of Bessel functions 
evaluated at the center, r=b, and free edge, r=a, of the circular plate. As a 
final step in the REDUCE program the driving point impedance was 
written in FORTRAN type statements to allow a simpler generation of a 
FORTRAN program. The REDUCE programs used in this study are in 
Appendix A. 

The flow chart for the numerical program used to calculate the 
driving point impedance is shown in Figure 13. The numerical language of 


FORTRAN was used in this study. 
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Calculate Young's Modulus 


Calculate Dimensionless 
Parameters 
Eq (2.22) and (2.23) 


Calculate Bessel Functions 


Calculate Driving 
Point Impedance 











Upper Frequency 
Limit Reached ? 






Figure 13. Flow Chart for FORTRAN Program. 


The FORTRAN program was written to so as to minimize the 
number of inputs required in order to calculate the driving point impedance 
over a range of frequencies. It was decided that the inputs would be the type 


of material, geometry, and temperature of the plate, and the upper and 
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lower frequency bounds of interest. The remaining information needed to 
calculate the driving point impedance would either be calculated or stored 
in table faction. 

As was noted by Irie, Yamada, and Aomura [Ref.12], the parameters 


592 and 532 can become negative such that the argument of the Bessel 


function is a complex value. In their study, to get around this, whenever 
522 or 53% become negative, the assumed solution went from the Bessel 
function of the first kind to the modified Bessel function of the first kind and 
the argument of the Bessel function was made into a real quantity. When 
dealing with viscoelastic plates, where the argument of the Bessel functions 
is always complex, this strategy was not acceptable. It was pointed out by 
McLachlan [Ref. 15, p. 27] that the Bessel functions of the first and second 
kind are still solutions to the Bessel's partial differential equation even with 
complex arguments. This means that the solution can be calculated if the 
Bessel functions are in complex notation. All of the Bessel function 
algorithms found on the Naval Postgraduate School's IBM 3033 are written 
for real valued arguments, not complex arguments, so could not be used. A 
computer algorithm had to be written to find the values of the Bessel 
function of the first kind, second kind, and their first and second 
derivatives. It was decided that the simplest method, and most straight 
forward, though not the most elegant, would be to write the algorithm with 
the Taylor series expansion of the functions. It must also be pointed out 
that this algorithm was written for the special case of Bessel function order, 
n, being zero. The Taylor series expansions used to calculate the Bessel 


functions and their first and second derivatives are [Ref 21]: 


3l 





>) Se (2.36a) 
92 9%. 42 9? 42. @? 2 2 cee 











3 5 7 
J Cx) = ee = Eee (2.36b) 
: 2°. 4 O°nde NiGun o aman heme 
2 4 6 
0:4 ee (2.36c) 
2:4 2-4-6 2:4 6 -8 
ez x 
¥ (x) = x In( a) +y hd o(®) 
(2.36d) 
il x4 0+9) oe 
Tt 9? 92. 4? 2 92. 42.6 ? 2° 3 
Y (x)= 2 In +y¥hT (x + $3 (0 
a : (2) J : (2.36e) 
aE X 1 X ( ‘leaietel | 
+33 = (| 1 + 2 aS eee 
1m) 9 3 3 97. 42. ¢ 2 _ 
4 2? xX 7 2 1 
Y 0) = 5] Inf ye) (ee ea 0x) 
0 hs (3) } 0 x" 0 =e (2.36f) 
2 4 
tea Gta) ee bata | 
2-4 2-4-6 
where x is a dummy variable and y is Euler's number which is equal to 
y= lim (1+5+5+ +4 -In(n)) (2.3) 
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For this study Euler's number was approximated by: 


y = 0. 5772156649015328606065. (2.38) 


The code used to calculate the Bessel functions was written in quadruple 
precision complex FORTRAN, which is an extension of FORTRAN 77 found 
on the Naval Postgraduate School's IBM 3033 computer. Using this 
precision of code for the Taylor series expansion insured that the resultant 
values had an accuracy of at least six significant places. 

In order to validate the programs written for this study, the driving 
point impedance of elastic plates would also be calculated and compared 
experimentally. This did not present any problems in the use of the 
programs, because an elastic material's Young’ modulus and shear 
modulus can be represented by complex values with zero imaginary part. 
The only addition required to the programs is the elastic material's 
characteristics. For the elastic material aluminum was chosen, and for 
the viscoelastic material LD-400 was chosen. 

It has been pointed out by Mindlin [Ref.10] ,and Reismann and 
Pawlik (Ref. 22] that the value of shear coefficient, k in equation (2.24), is not 
fixed but depends on the material's Poisson's ratio and the frequency of the 
first antisymmetric mode of thickness-shear vibration. Equation (2.24) is 
based upon only the first antisymmetric mode of thickness-shear vibration 
and corresponds to a Poisson's ratio of v= 0.176. For both the elastic 
material and viscoelastic material under study the Poisson's ratio is 


greater than this. It was decided that two values of the shear coefficient 


will be used for each type of material being studied. The values used are 
shown in Table 1. 

The material Young's modulus, or shear modulus, and the 
Poisson's ratio must also be supplied. For the aluminum plate this 
information was taken from Timoshenko [Ref. 23]. For the viscoelastic 
plates, the information was taken from Lee [Ref. 5] , which was based upon 
data from Nashif, Jones, and Henderson [Ref. 1, pp. 390-392], and from 
material data sheet supplied by the manufacture, United McGill [Ref. 24]. 
Only the data from United McGill (Ref. 24] was temperature compensated 
in the program written for this study. The information from United 
McGill [Ref. 24] was provided on a log-log scale chart which allowed 
determination of complex Young's modulus at three different 
temperatures, 68, 86, and 104 degrees Fahrenheit. For each temperature 
the complex Young's modulus was determined at four frequencies, 10, 100, 
1000, and 10,000 Hz. The information was taken in terms of the log of the 
value, so that the data was better conditioned for use in the numerical 
calculations. A four point, third order polynomial curve fit [Ref. 25] was 
applied at each of the above temperatures to allow calculating the modulus 
at other frequencies. At plate temperatures between those taken, a linear 
interpolation was used. The polynomial equations used are shown below. 


For a temperature of 68 degrees Fahrenheit: 


Die = 2.025 —Usdo2501. 7 0. 1000f7 = (. 01250f° 
E' log log log 


2 3 
Dama i 7s — OC O70S4t. 0.11 25f — 0. 01667f 
E" log log log 


For a temperature of 86 degrees Fahrenheit: 


2 3 
D_,=1.000+0.6125f, -0.1875f, +0.02500f 
E log log log 


E 


= 


For a temperature of 104 degrees Fahrenheit: 


D_ =0.9000 +0.08750f. +0.01250f7 +0.0000001590f° 
E' log log log 


2 3 
D_ =0.8000 +0.2542f -—0.06250f +0.008333f 
E" log log log 


where flog is related to the frequency by: 


f =log(f) 
log 


2 3 
D_, =1.050+0. 3875f, -—0.1000f, +0.01250f 
log log log 


(2.39a) 


(2.39b) 


(2.40a) 


(2.40b) 


(2.41a) 


(2.41b) 


(2.42) 


DF' and DE" are dummy variables which are related to the real part, E’, 


and the imaginary part, E", of the complex modulus by: 


(Pg) . | (2.43 a) 


D 
(PE) . | (2.43b) 


where E’ and E" are in the units of pounds per square inch. The material 
characteristics for the three reference temperatures are shown in Figures 
14 to16. The linear temperature interpolation was done in terms of DF' and 
Dk" and then converted into E' and E" by equation (2.43). It was felt that 
aluminum would be insensitive to temperature changes and the data from 
Nashif, Jones, and Henderson [Ref. 1, pp. 390-392] was too coarse to allow 
accurate calculations other than at room temperature, 75 degrees 
Fahrenheit. The shear modulus and loss factor equations are repeated 


here for the data from Nashif, Jones, and Henderson [Ref. 1, pp. 390-392]: 


G’ = 0. 00002503 -f °- 0.1752 -f °+ 457.5883 -f +29280 (244a) 


| F 1. sh 
-Q 52732 og( &) 


The moduli and loss factor for the viscoelastic material, from both 
sources are shown in Figures 17 to 19 for a temperature of 75 degrees 
Fahrenheit. This shows for the same material, there is little agreement on 


the material characteristics by different sources. 
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Figure 14. Viscoelastic Material Characteristics from 
United McGill [Ref. 24] at 68 Degrees Fahrenheit. 
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Figure 15. Viscoelastic Material Characteristics from 
United McGill [Ref. 24] at 86 Degrees Fahrenheit. 
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Figure 16. Viscoelastic Material Characteristics from 
United McGill [Ref. 24] at 104 Degrees Fahrenheit. 
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Figure 17. Comparison of Viscoelastic Material Young's Modulus 
from Nashif, Jones, Henderson [Ref. 1, pp. 390-392] and 
United McGill [Ref. 24]. 
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Comparison of Viscoelastic Material Shear Modulus from 
Nashif, Jones, Henderson [Ref. 1, pp. 390-392] and 
United McGill [Ref. 24]. 
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Figure 19. Comparison of Viscoelastic Material Loss Factor from 
Nashif, Jones, Henderson [Ref. 1, pp. 390-392] and 
United McGill [Ref. 24]. 


42 


Other inputs to the FORTRAN program are the density and the 
temperature of the plate. The summary of the constants used are listed in 
Table 1. A copy of the FORTRAN programs used can be found in 
Appendix B. The driving point impedance was calculated at 100 point 
across the frequency range of interest. For the case of the elastic plate, the 
theoretical impedance at a mode will tend to infinity. However, for the 
frequency range of interest, 100 to 2000 Hz, the chances of having the 
frequency of interest being the same frequency as a mode are very small, 
and an overflow of the computer is unlikely. The results, for all materials, 
are shown in the real and imaginary parts of the driving point impedance 
due to the waveguide predictions being based upon the real part of the 
impedance. In an effort to check the calibration of the experimental portion 
of this study, the driving point impedance was calculated in the range of 10 
to 105 Hz. The lower limit, 10 Hz, was set by the material characteristics of 
United McGill [Ref. 24] not going less than this, while the upper limit 
corresponded to the lower limit of the normal frequency range with a 5 Hz 


overlap. 
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TABLE 1. DIMENSIONS AND PROPERTIES OF 
MATERIALS UNDER STUDY. 


ELASTIC | VISCOELASTIC 
LD- 400 —DM-375 


Equations 
(2.39) to 





3.06713 W7iL101 3 170766 
x 10°32 x10 x 10-3 


The results for the elastic plates, shown in Figures 20 and 21, show 
that the plate is acting as a lumped mass in the low frequency range of 10 to 
105 Hz. A lumped mass has no real part to the impedance, and the 
imaginary increases at a linear rate, which is what is happening here. 
The results for the high frequency range, 100 to 2000 Hz, do not provide 
much information, except the frequency of the first mode, as shown in 
Figures 22 and 23. 

The results for the viscoelastic plates, shown in Figures 24 to 27, 
show that there is a significant difference in the driving point impedance 
due to the material characteristics used, either Nashif, Jones, and 
Henderson [Ref. 1, pp. 390-392] or United McGill [Ref. 24]. It also shows 
that at very low frequencies the viscoelastic material is behaving as a 
lumped mass, with little real part and a linearly increasing imaginary 
part. The effect of the temperature compensation of the material 
characteristics supplied by United McGill [Ref. 24] is demonstrated by the 
results shown in Figures 28 and 29. This shows that the impedance 
decreases and becomes smother with increases in temperature. 

In all cases the effect of a different shear coefficient, x, is 
insignificant and a shear coefficient of x*=n2/12 (k=0.9069) ail be used for 


the remainder of this study. 
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The dominant factors for the viscoelastic plates which effect the 
calculated driving point impedance, are the source of material 
characteristics and plate temperature for a fixed plate geometry, and the 
geometry of the plate. Figures 74 to 77 and Figures 96 to 99, in Appendix C, 
show the results for the five and four inch plates. 
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Figure 20. Effect of Shear Coefficient on the Theoretical Real Part of 
the Driving Point Impedance for a 6 in Radius Elastic 
Plate in the Frequency Range of 10 to 105 Hz. 


47 


105.0 


a 
——e 


44 K= a) () 69 
R= 07! 


ic 


PICAL 


LEGEND 


EORETICAI 


KORE! 


I 
| 


| 71 
Zl 





wD 





“> 02 0°0 


0°8 0°9 0 
(NI/OAS-AGT) JONVAAdNI 


Figure 21. Effect of Shear Coefficient on the Theoretical Imaginary 
Part of the Driving Point Impedance for a 6 in Radius 
Elastic Plate in the Frequency Range of 10 to 105 Hz. 


2000.0 


1610.0 


1620.0 


430.0 


1 


1240.0 


NCY (HZ) 


1050.0 


FREQUE 


4B0.0 670.0 860.0 


0.0 


au 





100.0 





0'O000c 0-000! 0°O 0'O000T- 0'000c- 
(NI/OAS—AIT) JONVAAdNWI 


Figure 22. Effect of Shear Coefficient on the Theoretical Real Part of 
the Driving Point Impedance for a 6 in Radius Elastic 
Plate in the Frequency Range of 100 to 2000 Hz. 
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Figure 23. Effect of Shear Coefficient on the Theoretical Imaginary 
Part of the Driving Point Impedance for a 6 in Radius 
Elastic Plate in the Frequency Range of 100 to 2000 Hz. 
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Figure 24. Effect of Shear Coefficient and Material Characteristics on 
the Theoretical Real Part of the Driving Point Impedance 
for a 6 in Radius Viscoelastic Plate in the Frequency 
Range of 10 to 105 Hz at a Temperature of 75.0 Deg. F. 
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Figure 25. Effect of Shear Coefficient and Material Characteristics on 
the Theoretical Imaginary Part of the Driving Point 
Impedance for a 6 in Radius Viscoelastic Plate in the 
Frequency Range of 10 to 105 Hz at a Temperature of 75.0 
Deg. F. 
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Figure 26. Effect of Shear Coefficient and Material Characteristics on 
the Theoretical Real Part of the Driving Point Impedance 
for a 6 in Radius Viscoelastic Plate in the Frequency 
Range of 100 to 2000 Hz at a Temperature of 74.5 Deg. F. 
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Figure 27. Effect of Shear Coefficient and Material Characteristics on 
the Theoretical Imaginary Part of the Driving Point 
Impedance for a 6 in Radius Viscoelastic Plate in the 
Frequency Range of 100 to 2000 Hz at a Temperature of 
74.5 Deg. F. 
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Figure 28. Effect of Shear Coefficient, Material Characteristics, and 
Temperature on the Theoretical Real Part of the 
Driving Point Impedance for a 6 in Radius Viscoelastic 
Plate in the Frequency Range of 100 to 2000 Hz. 
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Figure 29. Effect of Shear Coefficient, Material Characteristics, and 
Temperature on the Theoretical Imaginary Part of the 
Driving Point Impedance for a 6 in Radius Viscoelastic 
Plate in the Frequency Range of 100 to 2000 Hz. 


56 


HI. EXPERUMWENT 


A. EXPERIMENT SET UP 

A total of 6 plates were tested during this study. Three were made 
out of aluminum, an elastic material, and three were made out of LD-400 
material, a viscoelastic material. All were 3/8 of an inch thick with a plate 
from each material at a radius of four, five, and six inches as listed in 
Table 1. The viscoelastic plates where made from twelve inch square tiles 
that are 3/8 inch thick, this limited the maximum plate size to a six inch 
radius. The other two plates were chosen so that there would be a definite 
change in characteristics for a change in geometry and so that the ratio of 
thickness to radius did not become too large and deviate from thin plate 
theory. To determine the driving point impedance of the plates, the 


equipment and configuration of Figure 30 was used. 
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In Figure 30, WR is for Wilcoxon Research and HP for Hewlett Packard. 
The Wilcoxon Research F4/F7 shaker and Z7 measurement head are an 
integral unit that was mounted on a stand for the impedance 
measurements. The driving signal to the shaker is provided by the Hewlett 
Packard HP-3562A dynamic signal analyzer via the Wilcoxon Research 
PA7 power amplifier and N7/N9 matching network. The force and 
acceleration measurements of the plate were provided by the Wilcoxon 
Research Z7 impedance head mounting base, and then amplified by a 
Kistler 503 charge amplifier, for the force signal, or by a PCB 462A charge 
amplifier, for the acceleration signal, before going to the HP-3562A dynamic 
signal analyzer. The HP-3562A dynamic signal analyzer was operated in 
the swept sine mode for all of the tests to insure that the coherence of the 
data was at, or almost at, one. A minimum of 400 sample points were 
taken when the frequency range of interest was 100 to 2000 Hz, and a 
minimum of 160 sample points were taken when the frequency range of 
interest was 5 to 105 Hz. These were only a minimum as the analyzer was 
set to provide auto frequency resolution, or take more sample points, when 
the incoming data showed a large change from the previous data point. 
The analyzer was also set up to control the amount of power going to the 
shaker so that the incoming signals were of a quality that would keep the 
coherence at one. In the frequency range of less than 2000 Hz, the data was 
taken with the F7 piezoelectric vibration generator electrically secured and 
all data gathered with the F4 electromagnetic vibration generator. This 
was done because the F7 was producing considerable ring in the force and 


acceleration signals, which was degenerating the quality of the signal 


o9 


going to the analyzer. The F7 is normally used for the high frequency 
measurements and the F4 for the low frequency measurements [Ref.26]. A 
photograph of the four inch radius plate mounted on the shaker is shown in 


Figure 31. 





inch Plate Mounted on Shaker. 


Figure 31. Photograph of Four 
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The easiest method of determining the driving point impedance of the 
plate from the information, force and acceleration, being provided to the 
analyzer, is to start with the plate's frequency response and then 
mathematically convert this information into the driving point impedance 
by math functions built into the analyzer. The frequency response, H(f) , as 


measured by the analyzer, is: 


G 
SG) 3 (3.1) 


where GFA(f) is the cross spectrum between the force signal, F, and the 
acceleration signal, A, and GFF(f) is the power spectrum of the force signal 


[Ref. 27]. The transfer function, H(f), is the ratio of acceleration to force or: 


A 
H(f) == (3.2) 


The acceleration is related to the velocity by: 


Are ee (3.3) 


By combining equations (3.3) and (8.2) and rearranging the resultant 


equation the driving point impedance becomes: 
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j@ 


=F) 


(3.4) 


In the analyzer the frequency response, H(f), was first multiplied by (jw)-1 
and then the reciprocal of the resultant function was taken to be the driving 


point impedance, Z. 


B. PLATE MOUNTING DESIGN 

Some way of attaching the plate to the shaker must be designed and 
used. For use as a waveguide absorber, the plate will have a hole drilled at 
the center of the plate for a bolt or screw. This idea was used as the starting 
point for the design of the plate to shaker mount. A mount was 
manufactured from steel rod two inches long with 3/8 inch threads on one 
end, the shaker end, and 1/4 inch threads on the other end, the plate end. 
The plate was held between two 1/4 inch nuts and washers and a lock 


washer as shown in Figure 32. 






Nut 
Lock Washer 


Figure 32. Steel Plate to Shaker Mount. 


This mount proved to be unsatisfactory in operation. Its weight and size 
were totally against it. It resonated right in the middle, from about 700 to 
1200 Hz, of the frequency range of interest, 100 to 2000 Hz. The data taken, 
when using this mount, was useless. A new mounting system was 
designed, using aluminum, to minimize the weight and size of the mount 


and is shown in Figure 33. 





wi “eel Attach to Shaker 
8 


Figure 33. Aluminum Plate to Shaker Mount. 


This mount was extremely light weight and has almost no effect on the 
plate data. There was sufficient friction produced in the attachment to the 
shaker, that the plate and mount did not vibrate lose during the data runs, 
and a lock washer was not required. As a means of checking the 
effectiveness of this mount, the driving point impedance of the shaker and 
mount, with no plate, were determined over the same frequency range as 
with the plate. For these runs, the analyzer had the same configuration as 
with plates except that the auto frequency resolution was turned off. Ifthe 
auto frequency resolution were left on, the data runs would have taken too 
long, greater than 2 hours. The results of these runs showed that the 
shaker and mount were acting as if a pure mass of 1.48656 mslugs. The 


shaker and mount mass is much less than the mass of the smallest plate 


under study, so its effects are ignored. A 3/8 inch hole was drilled in the 
center of each plate to allow mounting. The plate sat on a 5/8 inch outside 
diameter boss on the shaker, and then was clamped down by the aluminum 
mount. Both the shaker and mount have an outside diameter of 5/8 inch 
which was used in all of the theoretical calculations for the value of inner 


radius, b. 


ce: VERIFICATION OF IMPEDANCE SCALE 

In all of the data runs, the calibrated output from the force and 
acceleration transducer, plus the gain of the charge amplifiers, was 
inputted into the analyzer to insure that the results of the analyzer were of 
the proper scale. As a means of checking the scale of these inputs, the 
driving point impedance of an elastic plate can be compared experimentally 
and theoretically at a certain frequency. At frequencies less than the first 
mode, the elastic plate can be considered as a lumped mass as shown in 


Figure 34. 


Displacement, x 





Mass, m 


Figure 34. Lumped Mass Model of Elastic Plate. 


Summing the forces on the mass, results in the equation of motion: 


F —-ma =F —-mx=0 


The acceleration, a, in terms of velocity, V, is: 


a= = joV 


The force, F, in terms of the velocity, V, 1s: 


F =jomV es 


The impedance, Z, then becomes: 


(3.5) 


(3.6) 


(3.7) 


(3.8) 


This allows calculating the impedance of the plate, based upon its weight, 


or mass, and the frequency of interest. A frequency of 100 Hz was chosen 


for this comparison on the elastic plates under study. The 100 Hz frequency 


is less than the lowest mode of the plates under study. The results of this 


comparison as shown in Table 2. Only the imaginary parts of-the 


impedance are looked at, as the model shows there is no real part, only an 


imaginary part. 
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TABLE 2. IMAGINARY PART OF IMPEDANCE FOR 
ELASTIC PLATES. 





ES 
i. Impedance, Imaginary Part, Ibf-sec/in 






The results of this comparison show that the analyzer and the associated 


equipment have the proper impedance scale. 


D. DATA RUNS 

Each of the six plates used in this study, were subjected to two 
different frequency ranges at room temperature. One data run was from 5 
to 105 Hz, for the purpose of checking the operation and calibration of the 
equipment, and the other from 100 to 2000 Hz, the frequency range of 
interest. As an added measure both of the six inch, elastic and viscoelastic, 


plates were run from 100 to 2000 Hz at an elevated temperature of about 100 


degrees Fahrenheit to check the effect of temperature on Young’ modulus. 
The surface temperature of the plate was measured with an Omega 871 
digital thermometer with a type K surface temperature probe. A number of 
points on the plate were measured before and after the data run. It was 
observed that the shaker was adding heat to the plates during the data run. 
This resulted in a temperature gradient of about five degrees Fahrenheit 
from the center of the plate to the outer edge. The temperature of the plate 
was taken as the average of all of the readings. 

On all data runs the coherence remained at one for every sample 
point taken. The repeatability of the data was tested on six inch radius 
elastic and viscoelastic plates by repeating the data runs after being 
removed from the shaker and then reinstalled, and by running the swept 
sine run in the decreasing mode instead of the increasing mode. There 
was almost no difference in the driving point impedance for these runs and 


the normal data run. It is concluded that the data is repeatable. 


E. EFFECT OF THE SHAKER ON THE IMPEDANCE 
MEASUREMENTS 
The shaker and the mount have mass which contributes to the data 
measurements. The effect of this mass can be taken into account by 


considering the model of the measurement system shown in Figure 35. 
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ZMeasured 





Figure 35. Impedance Measurement Model. 


The measured impedance for this model is: 


é Mount 2 Plate 


7 Measured =& Shaker Z (3:9) 


+Z 
Mount Plate 
If the impedance of the mount can be ignored, then the plate impedance 


becomes the measured impedance minus the shaker impedance or: 


é Plate é Measured “ Shaker (Sao) 
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This calculation can be accomplished in the analyzer very easily. In order 
to determine the impedance of the shaker and the mount, data runs were 
taken of the shaker and mount together, and with the shaker alone. The 
results of these data runs, Figures 36 to 39, show that the shaker and 
mount are effectively a lumped mass with little real impedance and a 
linearly increasing imaginary impedance with frequency. This occurred 
over both frequency ranges, 5 to 105 Hz and 100 to 2000 Hz, tested. The mass 
of the mount was measured to be 0.34946 mslugs, while the mass of the 
shaker, from m=Z/q@, is 1.0374 mslugs. The mass of the lightest plate is 
30.0123 mslugs, which is almost 100 times greater than the mount mass 
and 30 times greater than the shaker mass. As such the mass effect of the 
mount has been ignored and the plate impedance found by equation (3.10). 
It was noted that the data below about 15 Hz was very hard to accurately get 
without the mass of a plate being present. The data in this region has a fair 


amount of scatter. 
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Figure 36. Experimental Real Part of the Impedance of Shaker, With 
and Without Mount, in the Frequency Range of 5 to 
105 Hz. 
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Figure 37. Experimental Imaginary Part of the Impedance of Shaker, 
With and Without Mount, in the Frequency Range of 5 to 
105 Hz. 
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and Without Mount, in the Frequency Range of 100 to 


Figure 38. Experimental Real Part of the Impedance of Shaker, Wi 
2000 Hz. 
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Figure 39. Experimental Imaginary Part of the Impedance of Shaker, 
With and Without Mount, in the Frequency Range of 100 


to 2000 Hz. 
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In Figures 40 to 43 the difference between the measured and 
corrected plate impedance is shown. There is some small difference in the 
case of imaginary impedance and no difference for the real impedance, the 
shaker appears as a lumped mass. However the difference between the 
measured and corrected impedance is not considered to be significant to 
the conclusions of this study, and as such, the data will be in the 


uncorrected form for the rest of this study. 
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Figure 40. Experimental Real Part of the Driving Point Impedance of 
a 6 in Radius Viscoelastic Plate, With and Without Shaker 
and Mount Correction, in the Frequency Range of 5 to 105 
Hz at a Temperature of 75.0 Deg. F. 
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Figure 4L. Experimental Imaginary Part of the Driving Point 
Impedance of a 6 in Radius Viscoelastic Plate, With and 
Without Shaker and Mount Correction, in the Frequency 
Range of 5 to 105 Hz at a Temperature of 75.0 Deg. F. 
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Figure 42. Experimental Real Part of the Driving Point Impedance of 
a 6 in Radius Viscoelastic Plate, With and Without Shaker 
and Mount Correction, in the Frequency Range of 100 to 
2000 Hz at a Temperature of 74.5 Deg. F. 
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Figure 43. Experimental Imaginary Part of the Driving Point 
Impedance of a 6 in Radius Viscoelastic Plate, With and 
Without Shaker and Mount Correction, in the Frequency 
Range of 100 to 2000 Hz at a Temperature of 74.5 Deg. F. 
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F. RESULTS OF DATA RUNS 

For the elastic plates in the low frequency range of 5 to 105 Hz, shown 
in Figures 44 and 45, the plates appear to lumped mass as expected, with a 
small real part and a linearly increasing imaginary part with frequency. 
In the normal frequency range, 100 to 2000 Hz, the results, Figures 46 and 
47, for the elastic plates hold little information, only showing the location of 
the first mode, as was the case in the theoretical portion. There was almost 
no difference in the driving point impedance for the six inch elastic plate 
when data was taken at room temperature and when taken at an elevated 
temperature, as shown in Figures 48 and 49. This confirms that Young's 
and shear moduli for the aluminum plate need not be temperature 


compensated. 
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Figure 44. Experimental Real Part of the Driving Point Impedance of 
a 6 in Radius Elastic Plate in the Frequency Range of 5 to 


105 Hz. 
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Figure 45. Experimental Imaginary Part of the Driving Point 
Impedance of a 6 in Radius Elastic Plate in the Frequency 
Range of 5 to 105 Hz. 
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Figure 46. Experimental Real Part of the Driving Point Impedance of 
a 6 in Radius Elastic Plate in the Frequency Range of 100 


to 2000 Hz. 
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Figure 47. Experimental Imaginary Part of the Driving Point 


Impedance of a 6 in Radius Elastic Plate in the Frequency 
Range of 100 to 2000 Hz. 
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Figure 48. Effect of Temperature on the Experimental Real Part of 


the Driving Point Impedance for a 6 in Radius Elastic 
Plate in the Frequency Range of 100 to 2000 Hz. 
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Figure 49. Effect of Temperature on the Experimental Imaginary Part 


of the Driving Point Impedance for a 6 in Radius Elastic 
Plate in the Frequency Range of 100 to 2000 Hz. 
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The driving point impedance for the viscoelastic plate in the low 
frequency range of 5 to 105 Hz is shown in Figures 50 and 51. The results 
for the normal frequency rage of 100 to 2000 Hz are shown in Figures 52 and 
53. There is a major effect on the driving point impedance by the 
temperature of the plate, as was shown in the theoretical calculations. 
This is shown in Figures 54 and 55. 

Although the accuracy of the data is effected by such things as the 
calibration of equipment and the actual plate geometry, it is dominated by 
the material characteristics used. The material characteristics are a 
design specification provided by the vendor and the material may have a 
significant statistical variation in the material properties. The plate 
temperature also effects the material conditions for a fixed plate geometry 
and a different temperature measurement technique may provide a slight 


increase in the accuracy. 
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Figure 50. Experimental Real Part of the Driving Point Impedance 
for a 6 in Radius Viscoelastic Plate in the Frequency 
Range of 10 to 105 Hz at a Temperature of 75.0 Deg. F. 
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Figure 51. Experimental Imaginary Part of the Driving Point 
Impedance for a 6 in Radius Viscoelastic Plate in the 
Frequency Range of 10 to 105 Hz at a Temperature of 75.0 
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Figure 52. Experimental Real Part of the Driving Point Impedance 
for a 6 in Radius Viscoelastic Plate in the Frequency 
Range of 100 to 2000 Hz at a Temperature of 74.5 Deg. F. 
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Figure 53. Experimental Imaginary Part of the Driving Point 
Impedance for a 6 in Radius Viscoelastic Plate in the 
Frequency Range of 100 to 2000 Hz at a Temperature of 
74.5 Deg. F. 
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Figure 54. Effect of Temperature on the Experimental Real Part of 
the Driving Point Impedance for a 6 in Radius 
Viscoelastic Plate in the Frequency Range of 100 to 
2000 Hz. 
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Figure 55. Effect of Temperature on the Experimental Imaginary Part 
of the Driving Point Impedance for a 6 in Radius 
Viscoelastic Plate in the Frequency Range of 100 to 
2000 Hz. 
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IV. COMPARISON 


A. ELASTIC PLATES 

For the elastic plates in the frequency range of 10 to 105 Hz the 
agreement between the theoretical results and the experimental results is 
very good as shown by Figures 56 and 57. In the frequency range of 100 to 
2000 Hz the only thing that can be compared is the location, in frequency, of 
the impedance peak, as shown in Figures 58 and 59. For the three elastic 
plates studied the experimental peak occurred at a lower frequency, by 
about 100 Hz, than the theoretical peak. In the high frequency range, the 
theoretical calculations entirely missed the spike in the real portion of the 
driving point impedance, due to the frequency resolution used in the 
calculations. 

It is concluded that the ideas and programs used in this study have 


been validated to be correct by these results. 
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Figure 56. Theoretical and Experimental Comparison of the Real Part 
of the Driving Point Impedance for a 6 in Radius Elastic 


Plate in the Frequency Range of 5 to 105 Hz 


105.0 


FREQUENCY (IZ) 


ayes 
eS 
—_= 
<a 
~~ 

ie 
eee 
= 
—_ 


o™ 
— 
reer ey 
— 
ea 
= 
=! 
rae, 


20.0) 


rl 
4 
AL 


Lit 
4h4 
ad 

4 

4 





THEORI 
XPER 


4 
4 
¥ 


0 





Vga 
i 


0°8 079 0'+ 0-2 0°0 
(NI/OIAS-491) JONVGGd NI 

Figure 57. Theoretical and Experimental Comparison of the 
Imaginary Part of the Driving Point Impedance for a 6 in 
Radius Elastic Plate in the Frequency Range of 5 to 


105 Hz. 
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Figure 58. Theoretical and Experimental Comparison of the Real Part 


of the Driving Point Impedance for a 6 in Radius Elastic 
Plate in the Frequency Range of 100 to 2000 Hz. 
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Imaginary Part of the Driving Point Impedance for a 6 in 
Radius Elastic Plate in the Frequency Range of 100 to 


Figure 59. Theoretical and Experimental Comp 
2000 Hz. 


B. VISCOELASTIC PLATES 

In all cases there was almost no agreement between the 
experimental results and the theoretical results when the material 
characteristics of Nashif, Jones, and Henderson [Ref. 1, pp. 390-392] were 
used. But, as shown by Figures 60 to 63, there is agreement between the 
experimental results and the theoretical results when the material 
characteristics of United McGill [Ref. 24] are used. 

The remaining discussion deals with theoretical results with 
material characteristics supplied by United McGill [Ref.24]. In the low 
frequency range the experimental results would break or change at a lower 
frequency than the theoretical results would. This effect is similar to the 
experimental first mode for the elastic plates being lower than the 
theoretical. In the high frequency range the experimental real part of the 
driving point impedance was always lower than the theoretical, while the 
experimental imaginary part of the driving point impedance was always 
higher than the theoretical. It is as if the theoretical calculations are a few 
degrees in phase off from the experimental results. Figures 64 and 65 show 
the importance of including the temperature in the determination of the 
material characteristics. The theoretical and experimental results have 
the same shape but are shifted by a small amount. It is felt that this shift is 
due to the uncertainty in the material characteristics used in the 


theoretical calculations. 
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of the Driving Point Impedance for a 6 in Radius 


Viscoelastic Plate in 
a Temperature of 75.0 Deg. F. 


Figure 60. 
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Figure 61. Theoretical and Experimental Comparison of the 
Imaginary Part of the Driving Point Impedance for a 6 in 
Radius Viscoelastic Plate in the Frequency Range of 5 to 
105 Hz at a Temperature of 75.0 Deg. F. 
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Figure 62. Theoretical and Experimental Comparison of the Real Part 
of the Driving Point Impedance for a 6 in Radius 
Viscoelastic Plate in the Frequency Range of 100 to 2000 Hz 
at a Temperature of 74.5 Deg. F. 
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Figure 63. Theoretical and Experimental Comparison of the 
Imaginary Part of the Driving Point Impedance for a 6 in 
Radius Viscoelastic Plate in the Frequency Range of 100 to 
2000 Hz at a Temperature of 745 Deg. F. 
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Figure 64. Theoretical and Experimental Comparison of the Real Part 
of the Driving Point Impedance for a 6 in Radius 
Viscoelastic Plate in the Frequency Range of 100 to 2000 Hz 
at a Temperature of 100.5 Deg. F. 
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Figure 65. Theoretical and Experimental Comparison of the 


Imaginary Part of the Driving Point Impedance for a 6 1 
Radius Viscoelastic Plate in the Frequency Range of 100 


2000 Hz at a Temperature of 100.5 Deg. F. 
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C. CONCLUSIONS 

Ike The driving point impedance of a circular viscoelastic plate 
can be predicted by the programs developed by this study by knowing the 
plate's geometry and the material characteristics. The shape of the 
theoretical and experimental results is the same, but are shifted a small 
amount due to the material conditions used. 

ap The effect of differing shear coefficients on the driving point 
impedance is immaterial and can be ignored. 

S. The driving point impedance of a viscoelastic plate is controlled 
by the plate geometry and material conditions. The material conditions are 
greatly effected by the temperature of the plate. 

4, For viscoelastic plates, the magnitude of the real and 
imaginary part of the driving point impedance at a mode are about the 


same. 
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V. RECOMMENDATIONS 


Based upon the results of this study the following recommendations 


are made for studies continuing in this area. 


iG 


Investigate the effect of twisting of the waveguide absorber mount on 


the driving point impedance. 


. In the past study by Lee [Ref. 5] a constrained layer beam waveguide 


absorber was used with some success. A constrained layer circular 
plate can be used as a waveguide absorber and its characteristics need 


to be studied. 


. The effectiveness of the present viscoelastic circular waveguide 


absorbers on vibration control need to be studied and compared to a 
surface damping treatment and a point mass, all of the same weight. 
Investigate whether the circular plate waveguide absorbers can 
provide more energy dissipation per unit mass or volume of damping 


device than other forms of vibration control? 


This study has shown that plate geometry and plate temperature 
effect the driving point impedance. A sensitivity study should be 
undertaken to further investigate this dependance of the driving point 


impedance to changes in plate geometry and plate temperature. 
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5. Investigate the ability to design a circular waveguide absorber to 
remove a set amount of energy, at a point on the structure, at a certain 
frequency band. The design should be based upon the driving point 
impedance of the waveguide absorber and driving point impedance of 


the structure at the attachment point. 
6. Investigate the ability to design and use multiple waveguide 


absorbers to dampen over a wide frequency range or temperature 


range. 
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APPENDIX A. 
A. REDUCE PROGRAM TO SOLVE A FOUR BY FOUR MATRIX. 


te THIS PROGRAM IS DESIGNED TO SOLVE A 4 BY 4 MATRIX; 
“4 MATRIX METHOD; 

we PROGRAM TCAGA; 

OFF NATS 

OUT TT1HS$ 

LINEGENGINIC/2)); 


Dl: =DET MATCCB11,B22,855,899), (C2l,C22,6257624)), 
(C31,032,033,C034),0041,042,043,C0449); 
D2:=DET MATC(C11,C12,C13,C14),(B11,522, 83355439) 
(C€31,C32,C033,034),0C41,0462,043,C044); 
D3:=DET MATC((C11,C12,€13,€14), (O21, C22, €C2ameza 
(B11,B22,B33,B446),(C41,C042,C043,C0449); 
D4: =DET MATCCCl1,C12,C13,C14$), (C2, C227 C2. e237 
(C31,C32,€3535,C59), (Bil 8227555045, 
(C31,052,C033,C0354),(C41,042,043,C04649); 
SHUT TT1H$ 


BYE; 


pes 
b~4 
© 


B. REDUCE PROGRAM TO DETERMINE THE DRIVING POINT 
IMPEDANCE. 


Yu KK HHH HK HK HK HK RK KKK KKK KK HK HK HHH KKK HK HK HHH KKK KH KHKK KK AK KK KK KARKHKKKKR; 
A PROGRAM TA6C; 

i SPECIAL CASE OF N=0; 

OFF NATS$ 

OUT TTI1F$ 

7 THIS PROGRAM ATTEMPTS TO FIND THE SIX(6) FACTORS REQUIRED T0; 

vs DESCRIBE THE WAVE MOTION OF A PLATE WITH FREE EDGES AT R=A, AND; 
% A DRIVEN CLAMPED CONDITION AT THE CENTER R=B.; 

yA SPECIAL CASE OF N=0; 

7, LINE SIZE TO ALLOW FILE STORAGE; 

LINELENGTH (72); 

% ADDED DIFFERENTIATION RULES FOR THE BESSEL FUNCTIONS; 

meoH OF DELTAL; 

OPERATOR JN1,JNIDF1,JNIDF2,JNI1DF53; 

FOR ALL R LET DFCJNICR),RI=FDELTAIXJNIDFICRIZA; 

FOR ALL R LET DFCJNIDFICR),RI=DELTALXJNIDF2CR)ZA; 

FOR ALL R LET DFCJNIDF2COR),RI=DELTALXJNIDFSCR)ZA; 

% YH OF DELTAIL; 

OPERATOR YN1,YNLDF1,YNIDF2Z,YNIDFS; 

moR ALL R LET DFECYNICR), RI=DELIALXYNIDFICR)ZA; 
FOR ALL R LET DFCYNIDFICR),RI=DELTALXYNIDF2CR)IZA; 
FOR ALL R LET DFCYNIDF2(R),RI=RFDELTALXYNIDFSCRIZA; 
Zee vw OF DELTAZ2; 

OPERATOR JN2,JINZDF1,JN2DF2, JN2ZDF3; 

BUR ALL R LET DFCJNZ(R),R)=DELTAZXJNZDFICRIZA; 
FOR ALL R LET DFCJN2DFICR),RI=DELTAZXJN2DF2C(R)ZA; 

FOR ALL R LET DFCJN2DF2(R),RI=HDELTAZXIJN2ZDF3SCRIZA; 

% Yu OF DELTA2; 

OPERATOR YHN2,YN2DF1,YN2DF2,YN2DF3; 

FOR ALL R LET DFCYHN2CR),RI=DELTAZ2ZXKYHN2DFICR)ZA; 

Por ALL R LET DFC YN2DF1CR),RI=DELTAZXYN2DF2ZC(RIZA; 

FOR ALL R LET DFCYN2DF2(R),RI=FDELTAZKYN2DFSCRIZA; 

meee SN OF DELTAS; 

MmeRATOR JH3,JHSDFL,JN3SDF2Z,JNSDFS; 

mom ALL R LET DFC JHSCR),R)=EDELTASXIJINSDFILCR)ZA; 

mime ALL R LET DFCJNSDFICR),R)=DELTASXJNSDF2ZCRIZA; 

Burm ALL R LET DFC JNSDF2ZC(R),RI=DELTASXJNSDF3CR)ZA; 

fey OF DELTA3; 

MEPERATOR YHS,YH3SDF1,YN3DF2,YN3DF3; 

more ALE R LET DFCYHSCR),RI=DELTASXYNSDFICRIZA; 

BUR ALL R LET DFEFCYN3SDFICR),RIEFDELTASXYHSDFE2ZCR)ZA; 

moe ALL R LET DFCYNSDF2Z(R),R)=DELTASXYNIDF3CR)Z4A; 

fer ECIAL CASE; 

N:=0; 

% ASSUMED SOLUTION WITH THE BESSEL FUNCTION OF SECOND KIND; 


W1:= CALXJNICR) + AGXYNICR))XCOSCNXTHETA); 
W2:= CA2XJN2(0R) + ASXYNZ2CR))XCOSCNXTHETA) ; 
W3:= CASXIN3SCR) + ABXYNSCR)IXSINCNXTHETA) ; 


% POTENTIAL FUNCTIONS; 
KSIR:= (CSIGMAI-L)XDFCHI1,R) + CSIGMA2-1)*DFCW2,R) + DFCW3,THETAD/R; 


ICO 


KSIT:= CSIGMAlL-LIXDFCW1,THETA)/R + CSIGMA2Z-1)*¥DFCNZ2,THETAIZR - 
DFCW3,R); 

W:= Wl + W2; 

% MOMENTS AND SHEARS: 

DRC DFCKSIR,R) + CHUZRIXCKSIR # DECKSIIS THEA 

DACCKSIR + DFCKSIT, THETA) ISR + MUXDFCKSIR- RO 

MRT:= CD¥C1l-HUI/2)*CCDFCKSIR, THETA) - KSITIZR + DECKS Die 

KX¥X2XGXHXCKSIR + DFOCW,R)); 

Rex C~GXHK( COD URC THETA Ro. 

A BOUDRY CONDITIONS, 

ye R = A, MR=0,MRT=0, AND QR=0; 

MRBCA:= SUBCR=A,MR); 

MRTBCA:= SUBCR=A,MRT); 

SUBCR=A, QR); 

B AT THE CENTER, KSITR=0,KSI1=07 WwW =" 0; 

% WTO 2S THE DRIVING PORNT DISPEACETEM:: 


= 
7 
wou 


kD 
—d 
al 


KSIRBCB:= SUBCR=B,KSIR); 
KRSITSCB:— “SUBCR—8, Kol; 
WBCB:= SUBCR=B,H) - WTO; 
ue SOLVE FOR THE CONSTANTS Al TO Aé6; 
AAL:= © =~ BG4XC3GXC23XC12 + BGOXC34XC22*C13 + BGXC33XC24xCiZ 
BSOxXCSS3%XC229C14 - BSXCS2xXC244C)] 35 + B4GxXC 527025 ere 


+ C4GXB3*C234C12 - C4GXB3XC22XC135 = C44xCS3%E2 aie 
+ C4GXAC3535%C22*B1 + C44XCS32XB2XC13 = C4GXC32%C72 ea 
- CGS4B3*XC24XC12 + CG3XBaxCZ22ZxXC19 + C4SKCSGxp2 el 
= CG3XCS54XC22KB)] = C43XCS2ZXB2xC1G 4 C4S5*CS2xC2G- re 
+ CG2XB3*C24RC135 -— CG2XBSxCZ35XC14 = C42xXC5GxbZ ele 
+ CG2XC34XC2Z3XB1 + C42XC33KB2XC14 - CG2xXC53SxXCZ4xiD 


CCGG¥XC3S5XC2Z22XC11 -— CGGXC35*¥C21*C12 - CGGXC352*%C25xC1] +t C446XC32Z5Czie 
Cid + CGGKC31XC235XC12 - C4GXC31XC22*XC13 = C43XC34XC22e~C tee 
C3G4xXC2iXC1Z + CG3SXC352KXC2GXC1l1i -—- C43XC32*C21XC14 = CGSXC314C24-e Re 
+ CG35XC31¥C22*C14 + C42XCSGRXC23%C11 = C42*CS5GXE21KC13 — esc ese- 
C29*C11] + CG2*C35xXC21KC149 + C42xC3lxC2GxClSa = C42xCsixC23 et 
CG1=C34XC23XC12 + CG41XC349XC22XC13 + C41XC334C24XKC12 = CSlXCa eee 
CLG = CGIXCI324C24XC15 + CG1XE32*C25eC iG 
AA2Z:= CBG*¥C3GXC23XC1l1 —- BGXC34*XC21¥C13 - B4GXC33XC2GXC11 + BS 
ACSSXCZIXCIG. + BGXC31*C247C135 = BGXCS1 XO2540Cl4) [eee 
BS*C23*¥Cl1 + C4GXB3*C212C15 + C4GXCSSxECxCIT ~Veag- 
C354C21XB1] = C4G*Cl1xB2xC15 + 6CGGxC5 ll xe2Zr El eee be 
XC2GxKCI1 = CGSXBSxC2Z14C19 = C4SxXC3Gx%EB2-Cll + Cas eace 
C2le51 + C43xXC31XB2*C14 —- C43XCG1l4C2Gx Gl — Coheed 
¥ClS + CG1L¥S8SACC3xC14 + CGI ACS4XE25C13 ] CGI xACs Gee 
Bl = €414C33x%E24C14 + CGIxCS3xC24<B1)7 CC44%0s5s<eae- 
Cll = C46XC334C2Z21xC12 - CG4XC32*XC23KC11 + CG44XC32xC2Z1 NE 13 eee aa 
C31KC23%C1i2 - C44KXC31xC22KC13 - C43XC549XC22xCl1 + C45 Caa Gar. 
Cl2 + C43xC32*C24XC11] - CGS*xC32kC214C14G — C9354C31*C 25 cela eae 
C31XC22¥C14 + C462xXC54xXC25xXCl]1 -— C42xC3GxC2lx<Clsy G47 Cie C2 
Cll + C&2*C33KC214C14 + C42¥C31XC2GR%C13 - C42XKC51*E25<0l 4 Gee 
C34XC235*¥C12 + CG1*¥C54*C22¥C13 + C41XC33XC2GxCl2 = C4lxCSame re 
ClhG@. = C41 AC32xC26xCls + Cal es 2%Ccs eG), 
AAG:= € - BGXC34XC22xC1l1 + BGXC3G*C21*C12 + BGOXC32xC24*C11 
BGXC32XC21*XC14 - BG4XC31®C2GXC12 + BSXC31¥*C22XC14 


+ C4GXBS5*XC22xXC1l1 — CGGXBS3xXC2ZIXCI 2 = CoGxCicxE eee 

+ CGGKCS2KC2Z1AB1l + CGGXC31¥*¥B2*C12 - C44XC31*C22X%B1 

= C42XB3S4C24xXC11 + CG2XBSXC2Z214E1G6 + COZKCS4x82 Cr 

= C92KXC34GXC21*B1 - CGZXCSIAXABCXUIG tyCS2xCs14G24 ie 

+ CS1XB3XC24*4C12 - CG1XB3XC2Z2KXC14 - CG1XC34XB2XC1 2 

* CG1XCS54XC22XB1 + CGIXC32XB2XC1G = C6PRCS2xC2G hI 
(C4GXC33xXC22XC11 - CGGXC33*C21KC1l2 ~- CGGXC52*C23xXC11 + C4qxC32~e ee 


TZ 


feo t+ ©44GXC51XCZ25XC12 — C44*XC31XC22XC13 - C43XC34XC22XC11 + CG3xX 
C34*%C21*¥C12 + CG3XC32XCZ24G*KXC11 - C43XC32XC2Z1*XC1G - CG3XC31*XC24XC12 
+ CG35xXC31*¥C22*¥C1G + CG2*%C3GKC23XC11 - CG2*XC34*XC21*XC13 - CG2xXC33xX 
meGcxCl! + CG2xXC33xXC2Z1XC1S4 + CG2XCSIXCZ24GXC13 - C42XC31XC23*KC14 - 
CG1XC3GXC23XC1I2Z2 + CG1¥C34XC22XC13 + CG1¥C33XC2GXC12 -— CG1*¥C33XC22xX 
meg = CG1xXC52*XC24XC135 + C41XC3S2XC23XC14); 
BA5:= CBGXC3S5XC22XxCI]1 — BGXC3SXCZ1XC1l2 -— BGXCS2ZXCZ3XC11 + BG 

RCS2XCZIXC13 + BGXC3ILXCZ23XC12 - BGXC31xXC22xXC13 - CG3xX 

B3XC22XC11 + CG3X¥B3XC21XC12 + CG3¥C32*XB2XC1l1 - C43X 

C32*CZ1XBl1 —- C43XC3I1XB2XC12 + CG3XC31XC22XB1l + C4G2XB3 

mmesxCcl) = CG2*XBSXCZI1XC13 - CG2XC3SXB2XC11 + C42XC33xX 

Seep + CGZXC3IXBZ*XC13 - C42*¥C31*XC23XBl1l - CS1*XB3*C23 

meme + CG1XB3XC22XC13 + CG1XC33SXB2XC1l2 - CG1XC53*C22x 

me = CGIXC32XB2XC13 + CG1XC32XC2S5XB1)I/ 0 C44XC33XC22X 

Cll - CGG*XC33XC2Z21*XC12 -— CGGXC32XCZ3XC1l1 + C4GXC32XC2Z1XC13 + CG4x% 
MemreeCZSxCl2 —- CG6GXC31XC22*XC13 - CGSKXC3GXCZ2XCI11 + CG3XC34GXC21xX 
C12 + CG3XC32XC2GXC11 -— CG3KC32XCZ1XC14 - CG43XC31XC2GXC12 + CG3xX 
mec eexCl]G + CG2XC3SGXC235xCl11 - CG2XC3GXC2Z1XC13 - CG2*¥C33XC24X 
Cll + CG2XC33X*XC21*¥C14G + C42*C31*¥C2G*%C13 - CG2*%C31*XC23XC14 - CG1*% 
BoaaxC2sxCl2 + CG1XC34XCZ22XC13 + CG1XC33SXC2GXC1l2 - CG1X*XC33*XC22xX 
eg = CG1*XC352*¥C24XC13 + CG1XC32*XC235XC14G) ; 

feos = AA3; 

A6:= AAG; 

We DERIVATION OF Cll TO C44 AND Bl TO BG; 


Teo 


Cll:= SUBCA1=1,A2=0,AG=0,A5=0,MRBCA); 
C12:= SUBCA1=0,A2=1,AG=0,A5=0,MRBCA); 
C13:= SUBCA1=0,A2=0,AG=1,A5=0,MRBCA); 
C1G:= SUBCA1=0,A2=0,AG=0,A5=1,MRBCA); 
C21:= SUBCA1=1,A2=0,A4G=0,A5=0,QRBCA); 
C22:= SUBCA1=0,A2=1,AG=0,A5=0, QRBCA) ; 
C23:= SUBCA1=0,A2=0,AG=1,A5=0,QRBCA); 
C2G:= SUBCA1=0,A2=0,AG=0,A5=1,QRBCA); 
C31l:= SUBCA1=1,A2=90,A49=0,A5=0,KSIRBCB); 
C32:= SUBCA1=0,A2=1,A4G=0,A5=0,KSIRBCB) ; 
C33:= SUBCA1=0,A2=0,A4=1,45=0,KSIRBCB); 
C3G4:= SUBCA1=0,A2=0,AG=0,A5=1,KSIRBCB); 
CG1:= SUBCA1=1,A2=0,AG=0,A5=0,WT0=0,4KIBCB); 
CG2:= SUBCA1=0,A2=1,AG=0,A5=0,WTO=0,WBCB) ; 
CG3:= SUBCA1=0,A2=0,AG=1,A5=0,WTO=0,HBCB) ; 
CG4:= SUBCA1=0,A2=0,A4=0,A5=1,WT0=0,HBCB) ; 
Ba:= QO; 

Be: O; 

B3:= 0; 

B4:= WTO; 

Al:= AAI; 

A2:= [ya ae 

AG:= AAG; 

AD5S:= AAS; 


meet2:= INTCRXQR, THETA); 

moe-ll:= SUBCR=B,FOFT2); 

Beet: = SUBCTHETA=2XPI,FOFT1) -— SUBCTHETA=0,FOFTIL); 
Z:= FOFT/COMEGAXHTO); 

Pe orECIAL CASE OF N = QO; 

Of FORT; 

WRITE Z; 

OFF FORT; 

saul TTlF$ 

BE; 
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APPENDIX B. 
FORTRAN PROGRAM TO CALCULATE THE DRIVING POINT 
IMPEDANCE 


PROGRAM THESIS 
CK KK HHH HH HHH HK HK HK HH HHH HK HK HK HHH HH HHH HK HHH HK HHH HHH KK HH HH HK KH KKK HHH KK KK KKKX KKK 
Cx THE PURPOSE OF THIS PROGRAM IS TO GENERATE THE DATA HATCH. x 
Cx NECESSARY TO CALCULATE THE DRIVING POINT IMPEDANCE VS FREQUENCY x 
CH HH HK KK KH MH HEH KK HK HK HHH HHH HK HK HH HHH HH HK KH HH HH HK HK HH HK KKH KKH HK HHH KK KK KHKAKNX 


Cx LIST OF VARIA EES x 
Cx FL LOWER FREQUENCY OF INTEREST, HZ * 
Cx FH UPPER FREQUITEGY GF THERES lyn x 
Ce F FREQUENCY OF INTEREST. AZ ¥ 
Cx OMEGA ANGULAR FREQUENCY OF INTERES!, =2PIXF, RADYSEe ¥ 
Cx A RADIUS OF OUTER EDGE OF Finale. X 
Cx B RADIUS OF WAVE GUIDE ATTACHMENT, IN X 
Gx H THICKNESS OF VPEATE au X 
CX TEMP AVERAGE PLATE TEMP IRSDECGRESs-F x 
Cx IMA TYPE OF MATERIAL X 
Cx = 05> EEAS TIC x 
Cx = 1, VISCOELASTIC WITRSDATA FRUM- 6G Glee x 
Cx = 2, VISCOELASTIC WITH DATA FROM UNITED MCGILL X 
Cx RHO DENSITY OF PLATE, shuGe7 ix xs ¥ 
CH = YOUWG'*S MODULUS, Psl X 
Cx G SHEAR MODULUS, PSI X 
Cx MU POISSON'S RATIO X 
Cx PS f BESSEL FUNCTION OF THE FIRST KIND, EVALUATERM a X 
Gx DELTAXXY7A X 
Cx ) e.© BESSEL FUNCTION OF THE SECOND KIND, EVALUATEUMeT X 
Cx DEL TAAKYZA X 
Cz JXDELY FIRST DERIVATIVE OF THE BESSEL FUNCTION CF Re ® 
Gx FIRST KIND, EVALUATED AT DELTAXxY74a X 
ee YXADEIY FIRST DERIVATIVE OF THE BESSEL FUNCTION CF THe x 
Cx SECOHD KIND, EVALUATED Al Dela e172 X 
Cx 1 O5) aa & SECOND DERIVATIVE OF THE BESSEt FUNCTION Grete X 
Cx FIRST KID. EVACUATED AW Deni <7 A ¥ 
Cx YXDe2Y SECOND DERIVATIVE OF THE BESSED EUNCTIONSG Ess ¥ 
Cz SECOND KIND. EVAGUAEED AT DEE 12 X 
Cx DEUCES) DEMENSICHLESS PARANEDER * 
Cx DEGTAZ DEMENSIONHLESS PARAMETER X 
Cx Daas DEMEHSITUHLESS FARANEILER x 
Cx Seria. DEMENSIONCESS (PARSNERER X 
Cx SL Gitte DEMERSIUNLESS fora eer x 
e. DRIVING Pelil THPepAnce X 
Gx ZZ. COMPLEX TIME EDANGE ee br orc.) X 
Cx z MAGNITUDE OF ZZ, UBF-SEC7 ih x 
Cx ZR REAL PART OF ZZ, LBF-SEC7EN x 
Cex ZI IMAGINARY PART OF ZZ, LBF-SEC/7IN ¥ 
Cx ZDB Z AND X 
Cx K SHEAR COEFFICEINT X 
Cx Pal Sel Geer. X 
Cx x DUMMY ARGUMENT FOR BESSEL FUNCTION EVALUATION x 
CK RK HK KKK KK HHH KKK EK KKK EK KKK KKH KKK KKK KKK KA KKK KKK HHA HHHKKK HK KKK RX 


IMPELCIT. COMPDEXAXS “CAs D747 37)? 
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COMPLEXXS. -E5X%,6 

REAL MU,A,K 

CHARACTER ANS*%1, NAMEX6 

COMMONS BESSIAZJ1A,J1DF1IA,J1DF2A,Y1A,Y1IDFIA,Y1DF2A 
COMMON BESS1B/J1B,J1DF1B,J1DF2B,Y1B,Y1DF1B,Y1DF2B 
COMMON’ BESS2A7J2A,J2DF1A,J2DF2A,Y2A,Y2DF1A,Y2DF2A 
COMMON BESS2B/J2B, J2DFI1B,J2DF2B,Y2B,Y2DF1B,Y2DF2B 
COMMON’ BESS3A/J3A,JSDFIA,J3DF2A,Y3A,Y3SDFIA,Y3DF2A 
COMMON’ BESS3B/J3B,J3DF1B,J3DF2B,Y3B,Y3DF1B,Y3DF2B 
COMMON/GEOMTY/A,B,H,K 
COMMON/MAT/IMAT,E,G,MU,RHO, TEMP 
COMMON/CONST1/DELTAI, DELTA2, DELTA3, SIGMA1, SIGMA2 
PI= ARCOS(-1.0) 

HEITEC6, xX)" PI =*,P1 


INPUT OF MATERIAL TYPE 


5 CONTINUE 
Cree ERICMS C*CLRSCRN ') 
WRITECS,*)' THE FOLLOWING MATERIALS ARE SUPPORTED’ 
Hiker TECG sx) * : 


MAMAN 


VRREGG, x)" TYPE NAME CHARACTERISTICS FROM 
}! 

WRITEC6,%)! ’ 

Hike eco, x)" ELASTICCE) ALUMINUM GERE-TIMOSHENKO' 
WRITEC6,%)° : 

WRITEC6,%)* VISCOELASTICC(CV) LD-400 VIBRATION DAMPING' 
MRL LEC6G,x)' VISCOELASTIC(M) C-2204 UND RED MCGILL 


WRITEC6,%)' : 

Pee ECG, *)"' WHAT STYPE OF MATERIAL IS UNDER STUDY?! 

mel e(6,%)" TYPE™E, V, OR M.' 

READC5,*CA1)') ANS 

MIRITEC6,%)* ',ANS 

Pom Nee tt sGANUmeAtion NE. 'V"" SAND. ANS .NE. *M' ) GO TO 5 
IMAT = 0 

ieee ANS .EQ. 'V" ) IMAT 1 
Mec AHS .EQ. 'M* ) IMAT 2 
WRITEC6,*)" IMAT=°,IMAT 


mmo OF PLATE TEPERATURE 


6 WRITCC6,%)" WHAT IS THE AVERAGE TEMPERATURE OF THE PLATE IN DEGREE 
ieee CT EMP)" 
READ(5,*) TEMP 
peeve (G,*)' TEMP=', TEMP,' &* 
echt «Et eoo40 Oh TEMP ~-GT.- 104.0 ) GO TO 6 


otek @) 


BieUT OF PLATE GEOMETRY 


peeves O,™)" WHAT IS THE RADIUS OF THE PLATE IN INCHS?(A)' 
READ(5,%®) A 
WRITEC6,%)° A=",A,° IN! 
ee A .LE. 0.0) GO TO 5 
7 CONTINUE 


Olek @) 


C WRITEC6,%)" WHAT IS THE MOUNT RADIUS IN INCHS?(B)' 
C READ(C5,%) B 
fe 0.5125 


PrlneCG,*)' B=",B,* IN' 
Meee DLE. 0.0 ) GO TO 7 
8 CONTINUE 
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OO (06 f) Te OOO MOQ0O 9200) OOO (7 O)-6> 


2) 


oO OO 


O 


10 


15 


WRITEC6,%)* WHAT 15 THE THICKNESS OF THE PLATE INS ENGR Sram 


READ CS, <) 9H 

Hl = "5.07620 

WRITE(C6,%)* H=*n ee 
IF «CH .LE.~ O20 CGenoes 


INPUT OF FREQUENCY RANGE 


WRITEC6,%)*' WHAT IS THE LOW FREQUENCY OF INTEREST» IN HZ Gee 


READ(5,%*) FL 
WRITECO; 4%) Fla UR ees 
IF CFL 2LE. C20s > eGtr roar 


WRITEC6,%)°' WHAT IS THE HIGH FREQUENCY OF INTEREST, IN HZ? Geis 


READ(C5,%) FH 
HRITECG,*)° FH ea, ane 
Tr Cee CE. pie) OO ema 


INPUT OF LOCATION TO WRITE RESULTS 


WRITEC6,%*) * WHAT IS THE FILE NAME? ! 
READ(6,'CA)') NAME 


OPENCUNTT=15,FILE= NAME, FORM= ‘FORMATTED',ACCESS= "SEQUENTIAL’',ST 
LATUS= ‘NEW! 


TE OCCIMNAT “260 56>) =1 Orie = 299 
1F C IMAT .&@. 1 .ORS D4AT 2EQ. 2) 1POEn Tae 


START OF CALCULATION 
DO 90 I[=0,:iPOINT 
FREQUENCY OF INTEREST 


P= FL + CPs FlLIXREOA IC TZ ELOAT Chee) 
OMEGA = 2.0XPI<F 
WRITEC6,%)2" T=',1,' F=",F,*HZ OMEGA=", OMEGA, RaDZ see 


GENERATION CF DIMENSIONLESS PARAMETERS 
CALL CONSTCCMEGA, F) 
GENERATICGHN OF BESSEL FUCTIONS 


BESSEL FUNCTIONS FOR DELTA] AND RADIUS A 

X= DELTA 

WIRDVECGs.) * a 

CALL SESSELCXs 1A Ur 
PESSG. WFC PCS aU RD 

x > DELIA E77 A 

WRITEC6G,%4)' X=',X 

CALL BESSELCA JIB, JL DPI, JIDECB ey tea DOr ee 1 Gee 
BESSEL FURCTIONS FOR DELTAZ AND RAD] Joma 

Ks DEL ive 

WRITECGs <) “Kats x 

CALL BESSELCX,J2A,J2DF1A,J2DF2A,Y2A,Y2DF1A,YZ2DF2A) 
SOSEL FUNCTIONS FOR DELTA2 AND RADIUS B 

A = DELTAZ*BZA 

WRITEC6,%)* X=',X 

CALL BESSELCX, JZB, JZDF1B, J2ZDEZB, Y 252 2Ur bh eee 

BESSEL FUNCTIONS FOR DELTAS AND RADIUS A 

A == DELTAS 


t 


»JIDF2ZA,YIA,.YIDFIA, YIDF2ZA) 
TAL AND RADIUS 8B 


a 
m4 
1 
i= 


} 
fe 
a 
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WRITEC6,%)* X=',X 

CALL BESSELCX,J3A,J3DF1A, J3DF2A,Y3A,Y3DF1A,Y3DF2A) 
BESSEL FUNCTIONS FOR DELTA3 AND RADIUS B 

X = DELTASXB/A 

VIRITECG,*)" X=',X 

CALL BESSELCX,J3B,J3DF1B,J3DF2B,Y3B,Y3DF1B,Y3DF2B) 


IMPERANCE CALCULATIONS 


CALL IMPEDOCOMEGA,Z,2R,2ZI,ZDB} 
WRITEC15,20) F,2,2ZR,21,2DB 
90 CONTINUE 

WRITEC6,*)' DO YOU WANT TO CONTINUE?(Y/N).* 
READC5,'CA)') ANS 
Pelee Ri eMs CUCERSCRN I») 
Mees ANS .EQ. ty? )} GOTO 5 

20 eeu? Cemesernne 
t 
SUBROUTINE CONST COMEGA, F) 

CHE HEHE HK HK HEHEHE HE HEHE HEI HEHEHE HEE HE HE HEH HEE IE HEHE HEE HE HEHE HEH HE DEH HE I HEH HE HEHE HEH HE HE HE HE HE HE KH HH EH KH HH 


mOG@ (© © 


Cx THE PUPOSE OF THIS PROGRAM IS TO DETERMINE THE DIMENSIONLESS as 
CX PARAMETERS AND YOUNG'S AND SHEAR MODULUS FOR THE 
Cx VISCOELASTIC MATERIAL X 


CE HE EE EE HE KE HE EE EK EK EK EK HEE HEE HE IE EE EK EE EEE EE EE KE EE EE HK EE KK KKK X KR 


Cx LIST OF VARIABLE % 
Cx A PLATE RADIUS, IN “ 
Cx B WAVEGUIDE ATTACHMENT RADIUS, IN x 
Cx H PUATE THIECKHESS, IN ae 
Cx TEMP AVERAGE PLATE TEMP IN DEGREES F ¥ 
Cx IMAT tyPe UF SMA ERTLAL ¥ 
Cx =—0, ELASTIC MATERIAL ¥ 
Cx =i VESeCOELAITIC NITH DATA FROM G.G.LEE ¥ 
Cx = 2, VISCOELASTIC WITH DATA FROM UNITED MCGILL + 
Cx EF YOUNG*S MODULUS, PSI ¥ 
Cx G SHEAR MODULUS, PS! ¥ 
Cx ETHA Vis@UEPasollC EGss FACTOR B 4 
Cx GMAG VISCOELASTIC SHEAR MODULUS x 
CX D PeLKURAL RIGIDITY, PSI-IN x 
CX MU POTSS0N"’S RATIO ¥ 
Cx RHO Mass —UENSITyY GE PLATE, SLUGSZINX¥3 X 
eX K Sree cicoOerFIClI ENT ¥ 
CX OMEGA CIRCULAR FREQUENCY OF INTEREST, RAD/SEC ¥ 
Cx F FREQUENCY OF INTEREST, HZ ¥ 
Cx DELTAIL DEMENSIONLESS PARAMETER * 
ix DELTA2 DEMENSIGHLESS PARAMETER ¥ 
i DELTAS DEMENSIONMLESS PARAMETER % 
Cx SIGMAI DEMENSIONLESS PARAMETER % 
ex OAc DENENSIGHLEESS PARAMETER % 
ex DEETIZ DEDWTAl~x~2 % 
Cx DELT22 DES WAZx x2 ¥ 
CX DELT32 DELTA3**¥2 x 
Cx RR DEMENSIONLESS PARAMETER ¥ 
CX S DEMENSIONLESS PARAMETER ¥ 
Cx LS DEMENSIONLESS PARAMETER X 
Cx El ee be X 
CHE HE HE HK HE HE HEHE HE HEHE HEH HEHE EE HEE IEE IE HE HK HE HE IEE HE HE HE EH EIEIO EE EE EEE KKK 


COMPLEX¥8 DELTA1,DELTA2, DELTA3,SIGMA1,SIGMA2,L4,D,E,G,DELT12, 
Poe 22,DELT32 

REAL MU,K 

COMMON’GEOMTYZA, B,H,K 


ay 


COMMON/MAT/IMAT,E,G,MU,RHO, TEMP 
COMMON/CONSTI/ZDELTA], DELTA2, DELTA3,SIGMA1,SIGMA2 
PlsACost(=17 00 

HRITECX, X)* Pl = ee 


CALCULATION OF THE MATERIAL CHARACTERISTICS 
CALL YOUNGS(CF) 

CALCULATION OF FLEXURAL RIGIDITY 

Do = EXHX%S/ (12. 0X01 OU aMU 2) 

SHEAR COEFFICIENT 


K SQRTCPI*®*2712.0) 
K 0.9554] 

K 0.22700 

WRITECG,*)' K 
MRITECG, <7 2D 
Vinee, <) ' oo 
WRITEC6,%)' R 


DIMENSIONLESS PARAMETERS 


RR = CH/A)*®X%2/12.0 
WRITEC6,%)' RR=',RR 
S = (1.0706. 0X4Kxe 2x01, 0-MU) x CHAAD x x2 
HIRT TECGs 4) “3 5-"5 5 
LG4 = RHOXHXAKXGXOMEGAXXK2/(12.0%D) 
VIRITECG, 4)" LG=", 14 
DELTI12 = CLGXCRR+S#SQRT CCRR=530 4X2 4630714) 072.0) 
WRIT ECG ;x) ADE 2 elt 2 
DELT22 = CLGXCRRtS-SQRTC  RR-S)¥¥2+4.0/L9))72.0) 
Wind eCGs Xx) DEL co-7 so Dales 
DELTS2-= C2704 CRRMLG— 1 04507 (10D 
NRIRECG, <0 DELP S2= Debi 
SiGral-= DELT22 7 CKR¥U%—i3 075) 
RITE C6 4) * SiGMAl=', 5 8Giat 
SIGHA2 = DELTE2Z7(RRxUG=) 2075) 
WRITE(G, X)* SIGMAC=*,S1GMAZ 
Dowie = SGRTCDELTI 2» 
Pei eCG, <)* DELIA) ="s pebUAl 
Deli ae = oh) CDeel 22) 
HRELeCO, se) DElbwaZ-") DELTLAZ 
DELIAS = SORIC DED Ts2) 
Mel PEG, *)* DELTAS=", DELTAS 
G0 opie 0 
999 “CURIIHUE 
RETURK 
E:iD 
SUBROUTINE YOUNGSCF) 
OSS TESST SES SEL eS SPSS Se SesesSesesseSESSESL ESE S SESE ETS ESE SESE SESE SSE SES TS SS 


a 


CO) “Oo 0) Oe Om O10 Of QO G OO DOOD OOO O00 00) O00 Goo 


Cx THE MATERIAL CHARACTERISTICS X 
Cx THE PURPOSE OF THIS PRCGRAM IS TO GENERATE YOUNG'S MODULUS x 
Cx AND SHEAR MODULUS FOR THE MATERIAL UNDER STUDY x 
Cx C.D.HETTEMA G AUG 88 x 
CHK KK KH HK HK KK HH HHH HH HHH KEK KKK HK HHH HHH KH HK HK HIN HK KH HIKE KK HK KK EK EE KKK 
Cx LIST OF VARIABLES x 
Cx F FREQUENCY OF INTEREST; AZ X 
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PI Drak eens 


FLOG PoC rE? X 
T EMP pvenfoe PUATE TEMP IN DEGREES F X 
IMAT TYPE OF MATERIAL X 
= 0, ELASTIC X 

=Silveev FESCOELASTIC WITH DATA FROM G:.G. LEE X 

= 2, VISCOELASTIC WITH DATA FROM UNITED MCGILL ¥ 

RHO DENSITY OF PLATE, SLUGS/INXX3 x 
e YOUNG'S MODULUS, PSI ¥ 
G SHEAR MODULUS,PSI x 
ETHA LOSS FACTOR 
Gr REAL PART OF COMPLEX YOUNG'S MODULUS X 
EDP COMPLEX PART OF COMPLEX YOUNG'S MODULUS X 
EPX REAleraAnt OF COMPLEX YOUNG'S MODULUS AT TEMP X x 
EDPX COMPLEX PART OF COMPLEX YOUNG'S MODULUS AT TEMP X * 
DEP DUMMY VARIABLE FOR EP X 
DEDF DUMMY VARIABLE FOR EDP x 
DEPX DUMMY VARIABLE FOR EPX x 
DEDPX DUMMY VARIABLE FOR EDPX X 
MU POISSON *S RATIO X 
x 

ra 


CHK HK HHH KK HHH HK HK HK HK HH HK HK KKH KH HHH KKH HHH KH KH KK HHH KKH HH HK HK HK KKK KKK KKK KK KKK 


NINIAYMA ANY 


GOMPLEXXS £,G6 

REAL MU 

CHARACTER ANSX1 
COMMON/MAT/IMAT,E,G,MU,RHO, TEMP 


FOR USE IN UNITED MCGILL DATA 
FLOG = LOG1LOCF) 


KK KH KKK HHH HK KKK KKK KKK KKK KKH HK KKK KKH KK MK KKK HK KKH KKK HK KKK KHKKKKHK KKK KKKKKK 


ELASTIC MATERIAL DATA 


Meee IMAT .€Q. 0) THEN 
eee (11.05&6,0.0) 

MU = 0.33 

Gee E7(2.0%(1.0FMUD ) 
Reo! = 5.06/713E=-3 


CHK HK HK HH K HHH KH KH KH HK HK HK HK HHH HHI HK HH KH HEH HH HE HE HK HE HE HE HE HE HE HEH HE HE HE HEH EH HE HK KKK RK KKK 


NONNNY 


VISCOELASTIC MATERIAL CHARACTERISTICS FROM NASHIF, JONES, AND 
HENDERSON 


ebse TEC IMAT .EQ@. 1 2. THEN 
ee = OU OSXEXP(—-0 527 32X( CABS(LOGIOCF/62)) > xX*1.956)) 
GiaGe- 0 ,0000Z250S4F4X3 = 0.1752XFXX2 + G57.5883*F + 29280 
Gee  OnAGe((C].0,0.0) +.€0.2,1.0)XETHA) 

MU = 0.5 

eee GOxC2.0x( 1. 0+MU)) 

moe > 1. 7lIOlE-3 


CH HEH HK HE HHH HK HH HH HHH HHH KKK HHH HHH HH KKH MH HHH HHH IH HHH HH HH HHH HEH HK KAM K KKK KK KKK 


NNININNIH ANN 


VISCOELASTIC MATERIAL CHARACTERISTICS FROM UNITED MCGILL 
meoe IF € IMAT .EQ@. 2 ) THEN 
FOUR POINT, THIRD ORDER CURVE FITS FOR A SET TEMP 


fife OF 63 F 
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NMAGD 


MONADN 


OON 


OOM 


DEP. FOUR )=66F 
A3S==0-126999657e—0m 
Ac 0. 9999986 50 Ean 
Al=-U0..62500173¢6E— 07 
AO= 2.02999962 


DEPA = AO + ALXFLUG f AZXELUGXAZ + eAS PUG. 


DEDP FOR T=68F 
A3=-0.16666/7290E-01 
A2= 0.112500608 
Al=-0.7063551 P35E-on 
AO= 1.77500057 


DEDPA = AQ + ALXFLOG + A2XFLOGX*K2 + AZXFLOG*X3 


TEMP OF -SosF 


DEP FOR ls¢on 
AS= 0. 25000Z51S5E-01 
AdZ=—-O..TS 7501646 
Al= 0.612504184 
AQ= 0.999996543 


DEPB = AQ + ALXFLOG + A2XFLOGXX2 + A3SXFLOGXX3 


DEDESECK T=a6F 
AS- 021269997657 E701 
Ac=-0,.990999300cCE-O01 
Al= 0.387497604 
AO= 1.05000114 


DEDPB = AO + ALXFLOG + A2¥FLOG¥*2 + ASXFLOGRXS 


TEMP-0F- 1 0G: 16 


DEP FOR T=10GF 
AS= 0.158945682E-06 
oe. U1 cS tees boe=0 1 
Al= 0 ..6/50250076-01 
BU=" Uc 99 79 6aoo 
DEPC = AO + ALXFLOG 
DEDP FOR T=104F 


fo= 0.€33520616E-02 
A2=-0.6246990463E-901 
Al= 0.254164696 

AQ= 0.€00002154 

DEORC. <= -f0) eax LOG 


+ AZXFLOGXX2 + AISXFLOGXXS 


LINEAR PIT FOR SE? FREQUENCY 
TEM Cee ae 


le -C TEMP {Ge woe 70 


Also. 
DEP = <DEPS-DEra)*Cfeisr-con 007 le =) Eee 


A2¥FLOGEX2 + 


ASXFLOG¥X3 


8¢.0 ) THER 


DEDP =. COEDPB=DEDPAI~CTEMP=68 00716. 0 tee oee 


pe sage © ae Geet sti oe Gus Y 


DEP. = C(DEPC-DEr 3) CTE s=so 0 071c 0. epee 


E6.0 


AND. 


(EE 


sles 


104..0° 2) Then 


DEDP = CDEDPC-DEDPBIXCVEMNP-66.0)7 16 08 Jee 


ae Oe 


CONVERSION OF DATA TO PSI UNITS 


EP = 1.0E8X10**DEP¥1.45047E-5 
EDP = 1.0E8*1LOXXDEDP*1.G50G7E-S5 
FE = EPxXCi.0207 0) + EDPxCOl 0 iu 


PTA =" 2 DPZEP 
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MU = 0.5 
eee? Ce. OKC UU 
RaOms 1. /07S6Es5 
END IF 
RETURN 
END 
SUBROUTINE BESSELCXX, J0,JODF1,JODF2,Y0,YODF1,YODF2) 
CHE HE HE HE HK HE HE HE HEHE HEHE HEHE HE HE HE HE HEHE HEE HE HE HEE HE HE IE HK HEH HK HE HE HE DE HE HE HEE HEH I HE HE HEE EE EE KKK KK KKH 


Cx THIS PROGRAM CALCULATES THE BESSEL FUNCTION OF THE FIRST AND x 
Cx SECOND KIND ALONG WITH THE FIRST AND SECOND DERIVATIVE OF EACH X 
Cx fPaeo HIEL BE BY TAYLOR SERIES EXPANTION ¥ 
CX H = 0 x 
CE HE HE HE IE HE HE HE EE HE EE HE HE EE HE HE HE IE EE HE HE HE EE EE HE HEHE EE EK EK EE HE EE EHH EH EK KK KK HK KK KK RK KK 
Cx LIST OF VARIABLES X 
Cx M COUNTER 
Cx el COUNTER FOR DERIVATIVES, M1=M+1 x 
Cx 2M DOUBLE PRECISION M X 
Cx ZM1 DOUBLE PRECISION Ml x 
Cx ZA SUMBUR ST HESINVESSE OF Mf ¥ 
Cx x ARGUMENT OF THE FUNCTION ¥ 
Cx XXXD DOUBLE PRECISION VALUE OF Xxx * 
CX J0 BESSEL FUNCTION, FIRST KIND X 
Cx JODF1 Piro PeWERIVATIVE OF J0 ¥ 
CX JODF2 SECOND DERIVATIVE OF JO ¥ 
Cx JOA CHAHNIGE IN JO EACH STEP 
Cx JDFIA CHanGe IN JODFI EACH STEP ¥ 
Cx JDF2A CHANGE IN JODF2 EACH STEP X 
CX JOl ABS OF JOA x 
Cx JDEL1 ABS OF JDFIA * 
CX JDF2] ABS GF JDF2A X 
Cx YO BESSEL FUNCTION, SECOND KIND x 
Cx YODF1 Flrot DERIVATIVE OF YO ¥ 
Cx YODF2 SECONDS DERIVATIVE OF YO x 
Cx YOA CHANGE IN YO EACH STEP x 
Cx YDFIA CHANGE IN YODF1 EACH STEP X 
ex YDF2A CHANGE IN  yYODF2Z EACH STEP x 
Cx Y Ol ABS OF YOA ¥ 
Cx YDF11 ADowol 1 OELA x 
Cx YDF21 ABS OF YDF2A x 
Cx GAMMA EVEER’S CONSTANT X 
ex ERRCRM MAXIMUM ALLOWED ERROR X 
Cx ZZZMAX MAXIMUM SIZE NUMBER ALLOWED 
ex mel,C2Z,MC5 COUNTERS * 
CX D DUMMY VARIABLE 
CH HE HEE EK HE EEK EK HE HEH HE HE IEE HK EN I EK HE EK EH HK HE HE HE EH HEE HH HK HHH HK HEH KH KKK HK HK HHH KKK KK 


QUAD COMPLEX 
WoREICIT REAL¥16 CA-H,0-Z) 
DOUBLE COMPLEX 


tire lCl? COMPLEX*¥16 CA-H,0-Z) 

COMPLEX¥32 JOD, JODFID,JODF2D, JOAD, JDFIAD, JDF2AD,J01D,JDF11D,JDF21d 
COMPLEX%32 YOD,YODF1D,YODF2D,Y0AD,YDF1AD, YDF2AD,Y01D,YDF11D,YDF21D 
COMPLEX*32 X 

COMPLEX®3 JO,JODF1,JODF2 

COMPLEX*8 YO,YODF1,YODF2 

COMPLEXX*8 XX 

CHARACTER ANS] 


QNONMNMQ 
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OAD ANAYN 


MOON ANANYN 


60 


oxek@) 


OMY 


ANAND 


AOD 


Zen) = UDG 

ONE 1.0D0 

TWO ail) 

Pl = DARCOS( Oe) 

MAXIMUM ERROR ALLOWED 
ERRORM =) 170D=15 

MAXIMUM SIZE NUMBER ALLOWED 


ZZZMAX = 1.0D30 


GAMMA = 0.577215664901532860606512D0 


CONTINUE 
JO = ZERO 


CONVERSION OF X TO QUAD PERCISION 
xX = XX ® ONE 


CHECK TO PREVENT YO FROM BEING UNBOUNDED 


D= ABSCX) 
TP CD Saal 1 0D => 2 PE 


MIRITECG RO" XxAXR X 15 TOO SMAUI keer 


G60 70-59 
END IF 


CONTINUE 
JO CALCULATIONS 
Ml = M41 
= DBLECM) 
ZM1 = DBLECM1) 
SET UP OF TAYLOR SERIES 


De eee a. UO) 22 iets 


JOD = ONE 

JOAD = ZERO 

Ome = Ce 

JODF1D = -ONE#X/THO 
JDE IAD = 2eERU 
JODPLI = Ghee 
JODEF2D = -ONE/THO 
JDe2ZAS = ZERO 

JDe2 LD = OnEZ THOS 
Mono = 3 

MCT c= 1 

Cer =. 

MCS. =— 1 


EXPANSION OF TAYLOR SERIES 


ELSE DP uC Me Glewre: a ihe 
M2 = M72 


FOR J 


_ 


QAaAD 


ANN 


QNND 


ANN 


ANY 


20 


DA = ABS(JO01D) 

Lee Clem enkRORt ) THEN 

JOID = JOLDXCCX/ CTNOXZM) )*X2) 
JOAD = JO1D 

IF ¢ M2 .EQ@. MLAST ) JOAD = -J01D 


JOD = JOD + JOAD 
MGle= MC lod 
EhiD “LF 


FOR JODF 


DB = ABSCJDF11D) 

IF ¢€ D .GT. ERRORM ) THEN 

JDF11D = JDFI1LDXX**2/7 0 CTHOXCZM1-ONE) )¥C TWOXZM1 ) ) 
JDF1AD = - JDF11D 

reece .EQ. MLAST » JDFIAD = JDFIID 

JODF1D = JODF1D + JDFLAD 

Mc2 = MC2 + 1 

END IF 


FOR JODF2 
DC = ABSCJDF21D) 


fees D .GT. ERRORM ) THEN 
JDF21D = JDF21DXXXX2*(C TWOXZM1-ONED/CCTWOXCZM1-ONE) ) XC THOXZM1 


PO ACtHCxCZMI-ONE) =— ONE) ) 


JDF2AD = - JDF21D 

IF € M2 .EQ. MLAST ) JDF2AD = JDF21D 
JODF2D = JODF2D + JDF2AD 

feo = NCS + 1 

Si) LF 

MLAST = M2 

END IF 


CHECK IF SERIES CONVERGED 


D = ABS(JOD) 
DI = ABS(CJODF1D) 
= ABSCJODF2D) 
Bei Cl ow S00 mOpeme lr, Glew ZZZMAX .OR. DA .GT. Z2Z2Z2MAX .OR. 


molee.G). ZZZMAX OR. DB .GT. ZZZMAX .OR. 
Mig, Gi. ZEZMAX (OR. DE .Gi. ZZZMAX ) THEN 


MRITEC6,%)' *x*®x* FUNCTION HAS NOT CONVERGED *Xxx!' 
GO TO 50 

eo IF 

mee OA CLE. ERRORIT AND. M .GT. 0 .AND. 


moet. ERRORM SAND. DC .LE. ERRORM )} GO TO 20 


fa=- M+ I 

GO TO 60 

CONTINUE 

CONVERSION FROM QUAD TO SINGLE PERCISION 
JO = JODX1.0 

JODFl = JODF1D*1.0 

JODF2 = JODF2Dx1.0 

YO CALCULATIONS 


M = 0 
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oeke) 


QAOaAMY 


a 


2 G)G) 


AMON 


OOM 


AON 


SET “UP GPSTAY LOR @ sek res 
YOD = (CTHO/PI)DXCLOGCX/TWO)+GAMMA) XJOD 


YODFID = CTWO/PIDXCCLOGCX/ THOD+GAMMA)XJODFID + JOD/X) + X7PI 
YODF2D = (CTHO/PIDXCCLOGCX/TWO)+GAMMA) XJODF2D + TWOXJODFID/X - 


TJO07CK2X 2) PONE Z ES 
YOID = ONE 

TRE LD = x7 ini 
YDF21D = ONE/TWO 
on ERO 

MLAS 0 
mew 
eZ 
MES 


EXPANSION ROTATE GR SERIES 


Mea Mae ak 

Zh = DERE) 

ANTS Vii ot (OME ZZ 
M2 = hr 


FOR YO 


DA = ABSCYO1D) 

DEC De... Gl- ERRORM OS IREN 

YOID = YOLD¥CCX/CTHOXZM) 3 X¥2) 
OAD = =7 0 bp 

Le C2 bo. MEAS) 1 0A0e = 01D 


es 


i C1 IN 


YOD = YOD + THOXYOADXZMI/PI 

MGds aac tee 

EMO) i 

FOR YODF 

DB = ABSCyDF11D) 

ie DBY Oi. ERRORM Sone GTS 10 ine 
YDFIID = YDFLIDXXXx2/ CC TNOXCZM-ONE) XC TWOXZM) D 
YDELAD: = =YRELID 


iF CC 2 «60. MAST) oN DELAD S27 EEL» 
YODriD = YOCrID + 1DEIADA THe <2 Ars 
MG = bane 


FOR YODF2 


DC = ABSCrTDF2i0) 
Lf C DC «Gi wee RORG SAND il Gio 1 Oe tmieN 


YDF21D = YDF21DXXXX2xXCTWOXZM-ONEIZ CC TWOXCZM-ONE) ) XC TWOXZM) 
LAC THOR CZMSCHE = ONE)? 


YOR ZAD: = -=1DEZID 

PPC 42 26Q. (MLAST. DO YBR2ADG- Teal) 
YODF2ZD = YODF2D +°YDEZADXIWNOAZHIZE I 
MGiIe Sitter to) 


ENDS IE 
CHECK 2F SERIES CONVERGED 
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D = ABSCYOD) 

D1 ABSCYODF1D) 

D2 ABSCYODF2D) 

Ee Gaiewoles UO meOR sD eeGl ar ecZ2Znax  0R. DA GT. Z22ZMAX .OR. 

PD lereGiee 2 COMA Ke OR eb beso.  <ZZMAX . OR. 

ZUZeeor ICCC OK w Demerol. ZZZMAX ) THEN 
WIRITEC6,%)' %** FUNCTION HAS NOT CONVERGED ®xXxx! 
GO TO 40 
ti ae be 
LF CeDAS-LE. ERRORMS AND. M .GT. 0 ..AND. 

TDB .LE. ERRORM .AND2) DC LE. ERRORM >) GO TO 49 
GO TO 30 

40 CONTINUE 


CONVERSION FROM QUAD TO SINGLE PERCISION 


YO = YOD*1.0 
YODF1 = YODFID*1.9 
YODF2 = YODF2Dx1.0 
50 CONTINUE 
RETURN 
SUBROUTINE IMPEDOCOMEGA,2Z,2R,21,2DB) 
OSES CCST SES SSE SE SESE EES SS SESS SSS SSE SS SES ESE SESE SS STS SS SESS SS TESTS STS EES: 


olek@) 


CX Piece URE Ose Of ities sUBROM TINE TS]1O CALCULATE THE IMPEDANCE, x 
Cx AT THE DRIVING POINT, OF A CIRCULAR DISK BASED UPON BESSEL x 
CX FUNCTIONS, DISK GECMETRY, AND DRIVING FREQUENCY. xX 
Cx XXX SPECIAL CASE OF N=0 Xxx x 


CH HK HK HK HE HK HEH HEH EK HE HK EH EE HE EE KH KE HK EK HK KKH HK EH KEK KH HK HK HK KKK KKH KH KKK KKK 


CX EIST OF VARIABLE X 
CX ZL CUMRUEX IMPEDANCE, LBF-SEC/IN x 
Cx Z MAGNITUDE OF ZZ, LBF-SEC/IN 5 
Cx ZR Rel Pea Rien 225 —LBPSSeEC/ Ii % 
CX at IMAGINARY FART OF 22, LBF=SEC/7IN X 
Cx ZDB 2 IN DB, x 
CX A PLATE RADIUS, IN x 
Cx B WAVEGUIDE ATTACHMENT RADIUS, IN x 
CX H PLATE THICKNESS, IN * 
Cx IMAT TYPE OF MATERIAL x 
CX E YOUNGS MODULUS, PSI x 
Cx G SHEAR HORUIEUS:, Fol X 
Cx D FLEAURAEeRECIDITY, eSIl-1h ¥ 
Cx MU POISs0u oS RATIO x 
CX RHO MASSeBENSITY OCF PLATE, SLUGS/INXX3 x 
Cx K SHER COE Re Chen | xX 
Cz OMEGA CIRCUPS RO ERESUEICT OF INTEREST, RAD/SEC x 
eee @6dDELTAL DEE WS hOULESS PARAMETER x 
CX DELTA DENENHSIONLESS PARAMETER X 
Cx DELTAS DEMENSTONLESS PARAMETER x 
CX SIGMA] DEMENSIONLESS PARAMETER X 
CX SIGMA2 DEM SEGHEES> PARAMETER x 
O33 RR DEMENSIONLESS PARAMETER x 
Cx = DEMENSTUNEESS PARAMETER x 
CX L4 DEMENSIONLESS PARAMETER x 
CX tel Si: [Gees x 
CX ANSX INTERMEDIATE STEP X ANSWERS x 


CH HH HH HEH HE HEHE HE HE HEHE HE HE HE HEH HEH DED HE HE DE DEE HE DE HE HE HE HE HE DE HK DE DK HE HE DEH I HE HE HE EH HE HE EE HK HEH KH HH HH 
IMPLICIT COMPLEX*8 (€A,D,J,S,Y) 
COMPLEX*®8 2Z,E,G 
REAL MU,A,K 
COMMON’ BESSIA/J1A,J1DF1A,J1DF2A,Y1A,YIDF1A,Y1IDF2A 
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AQOAMY 


COMMON/BESS1B/J1B,JIDFIB,JIDF2B, 118; 7 TDP EE eee 
COMMON/ BESS2A/J2A,J2DF1A,J2DF2A,Y2A,Y2DF1A,Y2DF2A 
COMMON’ BESS2B/J2B,J2DF1B,J2DF2B,Y2B,Y2DF1B,Y2DF2B 
COMMON’ BESS3A/7J3A,J3DF1A,J3DF2A,Y3A,Y3DFIA,YIDF2A 
COMMON, BESS3B/J3B,J3DF1B,J3DF2B,Y3B,Y3DF1B,Y3DF2B 
COMMON/GEOMTY/A,B,H,K 
COMMON/MAT/IMAT,E,G,MU,RHO, TEMP 
COMMON/CONSTI/DELTA1,DELTA2, DELTA3,SIGMA1,SIGMA2 
PISACUS( =12.0) 


IMPEDANCE CALCULATION, FROM REDUCE PROGRAM 


ANSG=Y2DFIBXJ2ZDFIAXYIDFIAXJ1DFIB¥MUXSIGMA2**2- 
2. *¥Y2DFIBXJ2DFLAXY1DFILAXJI1DFIBXMUXSIGMA2* 
SIGMAL+Y2DFIBXJ2DFIAXYIDFIAXJIDFI1B¥XMU*¥SIGMAL 
¥¥2+Y2DFIBXJ2DFLAXYIDFIBXJIDFIAXSIGMA2*K2X 
SIGMALX¥DELTA1-Y2DFIBXJ2DF1AXY1DFIBXJ1IDF1AX 
SIGMA2X*X2XDELTAL-Y2DFIBXJZ2DFLAXYIDFIBXJIDFIAx 
SIGMA2XSIGMHAL¥*2x*DELTAL+Y2DFIBXJ2DFIAXY1DF1IBX 
JIDFIAXSIGMA2XSIGHALXDELTAL-Y2DF1BXJ2DFIA*Y1DF1B 
ASIDE LAXMUXSIGMAZCAX2+2 . XY CDFI BAS ZDE IAS IDELE 
XJLDFLAXMUXSIGMACXSIGMAI-Y2DFIBXJ2DFIAXY1DF1B 
¥JIDFLAXMUXSIGi4A1 ®¥2 


“ANS 3=Y2DFLAXJ2DF1B¥Y1DFELBXJ1DFILA¥SIGMA2*SIGMA1 


XX2*XDELTAL-Y2DFIAXJ2DFIB¥Y1LDFIBXJIDFLAXSIGMA2 
XSIGMALXDELTAL+Y2DFIAXJ2DFIBXYIDFIBXJ1DFLAXMU 
X¥SIGMA2**2-2 .XY2DFLAXJ2DFIBKYIDFIBXJ1 DFLAXMUX 
SIGMA2*¥SIGMA1L+Y2DFILA*XJ2DFIBXY1DFI1BXJ1DFIAXMU*X 
SIGMAL*®*2+Y2DFIBXJ2DFZA*Y1DFLAXJIDFIB¥SIGMA2 
¥X2¥SIGMALXDELTA2-Y2DFIBKJ2ZDF2AXYIDFIAXJIDF1B 
¥SIGMA2XSIGHALX*¥2XDEL TA2-Y2DFIBXJ2DF2ZAXY1LDFIAX 
JIDFIBXSIGHA2KXSIGMA1*DELTA2+Y2DFIBXJ2DF2AXYIDFIA 
XJIDFIBXSIGMA1L¥*¥2*XDELTA2-Y2DF1B¥J2DF2AXY1DFIB 
¥JIDFIAXSIGMAZKX2*SIGMALXDELTAZ+Y2DFIB¥J2DF2A 
XYIDFIBXJIDFIAXSIGHA2KSIGMALXX2XDELTA2+Y2DF1 BX 
J2DFZAXYIDFIBXJIDFIAXSIGMAZXSIGMALXDELTA2-Y2DF1B 
RIJ2ZDFZAXY1IDFIB¥JLDFIAXSIGItAL*®X2*XDELTAZ2-Y2DF1B 
XJ2DFLAXYIDF2AXJIDFIBKSIGMACKKOXSIGMALXDELTAL 
+Y¥Y2DFIBXJ2DFIAXYIDF2AKJIDFIBXSIGMAZKX2XDEL TAL 
+Y2DFIBXJ2DFLA¥YIDF2AXJIDFIBXSIGMALKSIGMA1 €X2 
¥DELTAI-Y2DFIBXJ2DFIAXY1 DF2AXJ1DF18*SIGMA2X 
SIGMAI*DELTA1L+ANS4 


“ANS2=-Y2DF2A*J2DFIBXY1LDELAXJ1DF1 B¥SIGMA2*X2X 


SIGMALXDELTAZ2+Y2DFZAXJ2DF1IBXYLDFIA¥ILDFEIBX 
SiGHASXSIGHAIXXZXDELTAZ+Y2CDEZAXJZOEIT xy Gee 
JIDFILBXASIGHA2*SIGHA]XDELTA2-Y2DEZAxXI2ZDEIEx pews 
RJIDFIBASIGIMALX¥2XDELTA2+Y2DFZAXJ2DFIBXYIDFIB 
XJIDFIAXSIGHAZ*¥2*S IGMALXDELTA2-Y2DF2AxJ2DF1B 
XYIDFIBXJIDFIAXSIGMA2*SIGHAIX¥2XDELTA2Z-Y2DF2AX% 
JZDFLEAXYIDOFIBXJIDFIAXSIGHNAZ*SIGMAL<BEt Ace yeni, 
XJ2DFLIBXVYIDFIBKXJIDFLAXSIGHAL¥R2XDELTA2+Y2DF1A 
XJZDFIBKYLDF2ZAXJIDFIBXSIGHAZXX2XSIGMALXDELTAL 
-~Y2DFlAAJ2ZDFIBXY1DF2ZAXJ1DFIBXSIGHAZK=2 Denia 
~Y2DFLAXJ2DFIBXY1LDF2AXJIDFIBXSIGMAZXS IGMA] X¥2 
X¥DELTAL+Y2DFIAXJ2DFIBXYIDF2AXJIDFIBXSIGIHA2* 
SIGHALXDELTAIL-Y2DFIAXJ2DFIBXYIDFIAXJIDFIBXMUX 
SIGHA2FK2+2. XY2DFLAXJ2DFIBXYIDFIAXJIDFIBXMUX 
SIGMAZ*XSIGMA1L-Y2DFLAXJ2DFIB¥XY1IDFIAXJ1DFIBXMUx 
SIGMAL¥X*2-Y2DFIAXJ2DFIB¥YIDFIBXJIDFIAXSIGHA? 
*¥¥2XSIGMALXDELTAL+Y2DFIAXJ2DFIBXY1IDFIBXJIDFIA 
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. XSIGMAZ*®X2XDELTAL+ANS3 

ANS1=2.XBXGXHXKX*¥2XPI XDELTAZXDELTA1LX¥ANS2 

ANSIL3=Y2BXJZ2DFIBXYLDFLAXJIDFILAXSIGMA2*SIGMA] Xx 
Z2X¥DELTAL¥¥2-Y2BXJZ2DFIBXYIDFLAXJIDFIAXSIGMA2* 
SIGMALXDELTA1L¥X2-Y2BXJ2DFIBXY1DFIAXJIDFI1AX 
SIGMA1LX¥X2*¥DELTALXX2+Y2BXJZDFIBXYIDFIAXJIDF1AX 

. SIGMALXDELTA1 ¥¥2 

Beene -2.XY2BXJ2DFLAXY1DFZAXJIDFIBXSIGMA2*X 

SIGMALXDELTA1LX¥¥2+Y2BXJ2ZDFIAXYI1DF2ZAXJIDFIBX 

SI GMAZ *UGeWAl*<2-1 20%) 20rlax1)] DEFLIAXJIDFIBXMUX 

SIGMAZHSIGMAISDEETAI+Y ZBEXIZDFIAXYIDFIAXJIDFI1B 

X¥MUXSIGMAZ2¥DELTAIL+Y2EBXJ2DFLAXYIDFILAXJIDFIBXMU 

XSIGMAL¥X¥2XDELTAL-Y2BXJ2DFIAXYIDFIAXJIDFIBXMU 

XSIGMALXDELTAL-Y2BX¥J2DFIAXYI1DFIBXJIDFI1Ax 

SIGMAZ*¥SIGMALX¥X2*XDELTA1L¥X¥2+2.*Y2BXJZ2DFIAXYIDFIBxX 

JIDELAXS IGMAZ*SUGMA TXDELTAIXX2-Y 2BXJZDFIAXYIDFIB 

TIER PAX ST GHAZCSDELTAISxX2Z+Y2bxJ2DFIAXYIDFIB 

. KSJLDFIAXMUXSIGMAZXSIGMAILXDELTAlL-Y2B#J2DF1AxX 

. YIDFIBXJILDFIAXMUXSIGMA2XDELTAI-Y2BXJ2DF1Ax 

. YIDFIBXJIDFIAXMUXSIGMA1L¥*2xXDELTAIL+Y2BXJ2DFIAX 
YIDFIBXJIDFIAXMUXSIGMAILXDELTAIL-~Y2BXJ2DF1 Bx 

YIDFZAXJIDFIAXSIGMAZ2XSIGMA1LXX2*XDELTALX*2+Y2BX 

JZOFIEXYIDEZAXJIDEFIAXSIGMAZXSIGMAIXDELTA1¥X2+Y2B 

XJZDFIBXYIDFZAXJI1DFLAXSIGMA1LX¥*¥2*DELTA1*X2-Y2B 

ICU Ey) vee. oP lax Shona] xDEL TAIXxXZ*ANS13 
“ANSIL1=-Y2DF1BXJ2BXY 1 DF2AXJ1 DFLAXSIGMA2XSIGMALX 

. DELTALXX2-Y2DFIBX¥JZ2BXY1DF2AXJIDFIAXSIGHA1XX2* 
DELTALX¥2+Y2DFIBXJ2BXY1DFZAXJIDFIAXSIGMAI1* 
DELTALX¥2-Y2DFIBXJ2BXY1DFIAXJIDFLAXSIGMA2xX 
SLCMALA* 2* Deby x2+Y2DelexJI2BXY I DELAXJIDFIAx 
SGine st CMa Copeman cect ec OP) BX IZEXYIDPIAXJIDFIA 
ASG MA lee XK DEP wA lee VCD ET BX ICE XYIDFIAXJIDFIA 
XSIGMALXDELTAL¥X2-Y2BXJ2DF2AXYIDFIAXJIDFIB 
X¥SIGMACXSIGMALX¥2*¥DELTAZXDELTAIL+Y2BXJ2DF2AXYIDFIA 
AJP IBS IT GMAZXStGMal*XDELTAZC*DELTAIL+Y2BxXJ2DF2AxX 
YIDFIAXJIDF1BXSIGMA]®X2*DELTA2Z*DELTAL-Y2BXJ2DF2A 
YP DPIAX JI DEB <srGriaIxXDEL TAZ*DEL TAI+Y2B*x 
V2DEZAXY IDFIBXILDPIAXSIGHAZXSIGMAIXX2*DEL TA2X 
DELIAL=Y2BxJ2DFCAXYIDFIBXJ] DFLAXSIGMAZXSIGHA] 
XDELTA2ZX¥DELTAL-Y2BXJ2DFZ2ZAXYIDFIBXJ1DFIAX 
SIGMAL¥X¥2XDELTAZXDELTAIL+Y2BXJ2DF2AXYIDFIBXJIDFIA 
¥SIGMAILXDELTA2XDELTAL+Y2BXJ2DFIAXYIDFZAxX 
JIDFIBXSIGHMAZXSIGMAlL¥X¥2*XDELTAL¥X2+ANS12 

“ANS10=~ Y2ZDrP lp J2OrPTAxy UDPZAxJ1TB*SIGMAZ=X2x 
SIGMA. =DELTAZHDELTAI+TYZDFIBSIZDFIAXY1DEZAXJ1B 
¥SIGHAZXAZ2XDELTAZXDELTAI+YZ2DF1B8XJ2DFIAXY1 DF2ZAx 
JIBXSIGHAZXSIGHALXDELTAZ2*DELTAI-Y2DFI1BXJ2DF1LAxX 
MODE CAxX IL EAS Grace Deb TACXDELTAL+Y2DFIBXJZDFIA 
XYLDFIAXJIS*¥MUXSIGMAZXX2XDELTAZ-Y2DFIBXJ2DFIA 

XY 1 DFIAXJILBXMUXSIGMA2XSIGMAILXDEL TA2-Y2DF1Bx 
J2DFIAXY1LDFIAXJIBXMUXSIGMAZ2*XDELTA2+Y2DF1 Bx 
JZDFLAXYIDFIAXJIBXMUXSIGMAI1XDELTA2+Y2DF1BxX 
J2DFLIAXYIBXJIDFIAXSIGMA2**2*SIGMALXDELTA2X*DELTAI 
mrelrlBxXJIZDFIAXY IBXJTBDFIAXSIGMAZCKXZ*¥DELTA2ZX 

DEE TA reUPILBXJ2DEVAXY IBXJIDEIAXSIGMHACXSIGMAL 
ADELTAZXDELTAI+Y2DFIBXJZ2DFIAXYIBXJ1DF1Ax 
SIGMAZ2XDELTAZXDELTAI-Y2DFIBXJ2ZDFIAXYIBXJIDFIA 
RiqUxol OMACKACH PEL IAetY 2ORPILB*JUZDFIAXYIBXJIDFIA 
AmMUXS CMA ZEST GMAIXDELTACTYZDFIBXJZDFIAXY1 BX 
JIDFLIAXMUXSIGMA2*DELTA2-Y2DFIBXJ2DFIAXY1I BX 


Pz 


JIDFIAXMUXSIGMALXDELTA2+Y2DFIBXJ2BXY1DF2AX 

» JIDFIAXSIGHA2XSIGMA1L¥*¥2*XDELTALX¥2+ANS11 

ANSO=Y2DFIAXJZ2BXY1DFLA¥ JI DFIBXMUXSIGMA2XSIGMA] 

» ¥DELTAI-Y2DFILAXJ2BXY1IDFLAXJIDFIBKMUXSIGMA2* 
DELTAIL-Y2DF1AXJ23*%Y1DFIAXJIDFIBXMUXSIGMA1Xx2x 
DELTAI+Y2DFIAXJ2BXYIDFIAXJIDFI1BXMUXSIGMAL* 
DELTAL+Y2DFIAXJ2BXYIDFIBXJIDFIAXSIGMA2XSIGMA] 
RX2XDELTALAX2-2.*XY2DFIAXJ2BXY1 DF1BXJ1DFIAX 
SIGMA2*SIGMAL¥DELTALXX2+Y2DFIAXJ2BXYIDFIBXJIDFIA 
XSIGHMACKDELTALXAZ=Y2ZDFIAXI2BXY IDPS ee 
XMUXSIGHAZKSIGMALXDELTAL+Y2DFIAXJ2BXY1DF1 Bx 
JIDFIAXMUXSIGMA2XDELTAL+Y2DFlAXJ2BXY1DF1 Bx 
JIDFIAXMUXSIGMAIXX2XDELTAL-Y2DFIAXJ2BKXYIDFIBx 
JIDFLAXMUXSIGMAL¥DELTAL+Y2DFIBXJ2DF2AXYIDFIAX 
JIBXSIGMACKX2XSIGMAIXDELTAZXX2-2 .XY2DFIBXJ2DF2AX 
YIDFLAXJIBXSIGMA2XSIGHMALXDELTA2*X2+Y2DFIBXJ2DF2A 
AYIDFIAXJIBXSIGMALXDELTA2Z*X2-Y2DFIBXJ2DF2A 
XYLBXJIDFIAXSIGMA2**2XSIGMAILXDELTA2%X2+2.x*Y2DF1B 
XJ2ZDFZAXYIBAJIDFLAXSIGMA2XSIGiHALXDELTA2%X2- 

Mae teatime a aaa) AT a 
S10 

“ANSS =“Y2DE1AxJ2DF1 B¥Y1DE2AxJ1B¥S1GMA2XSIGNAL* 
DELTA2*DELTAL+Y2DFILAXJ2DFIBXYIDF2AXJIBXSIGMA2 
*DELTAZ¥DELTAI=Y2ZDFIAXJ2DEI Bx; IDE Ax pie 
SIGMA2*XX2XDELTAZC+Y2DFLIAXJ2DFIBXY1DFIAXJ1BXMUX 
SIGMA2RSIGMALXDELTA2+Y2DFIAXJ2DFIBXYIDFIAXJ1B 
XMUXSIGMACKDELTA2-Y2DFIAXJ2DFIBXY1DFIAXJIBXMU 
¥SIGMALXDELTA2Z-Y2DFIAXJ2DFIBXY1IBXJIDFIAX 
SIGHACK*2RSIGMALXDELTA2XDELTAIL+Y2DFIAXJ2DFIBXY1B 
XJIDFIAXSIGMAZKX2XDELTA2XDELTAI+Y2DFIAXJ2DF1 Bx 
YIBXJIDFLAXSIGMA2XSIGMALXDELTA2*XDELTAI-Y2DF1Ax 
JZDFIBKXYIBXJIDFIAASIGMA2ZXDELTAZXDELTAI+Y2ZDEIA 
XJ2DFIBXYLBRJIDFIA¥XMUXSIGMA2*XX2xXDELTA2-Y2DF1A 
XJ2ODFIBXYIBAJIDFIAXMUXSIGMA2XSIGMAILXDELTA2- 
Y2DFIAXJ2DFIBXY1BXJI DFLAXMUXSIGMA2¥DELTA2+ 
Y2DFLAXJ2ZDFIBXY1BXJIDFIAXMUXSIGMAI1XDEL TA2- 
Y2EDFLAXJ2BXYIDF2AXJIDFIBXSIGHA2XSIGMALXX2% 
DELTALXX2+2.*¥Y2DFIAXJ2BXY1DF2AXJIDFIBXSIGMA2* 
SIGMAI¥DELTAIL¥X2-Y2DFIAXJ2BXY1DF2AXJ1IDF1Bx 
SIGHAZXDELTALXX2+AnS9 
“ANS7 =-Y2DF2A¥J2BXY1 DFLAXJ1 DF1B¥SIGMALX¥2X 

- DELTACXDELTAI+Y2DF2AXJ2Bx*Y1DFIAXJIDFIS*SIGMAL 
XDELTAZADELTAL-Y2DFZA*IZBXY LUPE UB =e PAs 
SIGHMAC#SIGMAILX*2XDELTA2ZXDELTAL+Y2DF2AXJI2BXY1DF15 
XJIDFIA¥SIGMA2XSIGMAILXBDELTA2XDELTA1+Y2DF2AXJ2B*¥ 
YIBFIBDBXIJIDFIA¥XSIGHAILXX2xDELTA2¥DELTA1L-Y2DF2AXJS2B 
AVYIDFIBXJIDFIAXSIGMALZXDELTA2*DELTA1L-Y2DF1AxX 
JEDFZAXYIDFIBMII BASIGUIAZAz 2x51 GitAl DECI ee ae 
YODFIAXJ2DFZA*YIDEIBXJIBXSIGHA2RX24DEL TALE 2+ 
Y2ZDFIA*J2DF2AXYIDF1BxIJ1exS1GighexSt Gia ~DeELi AZ 
¥X2-Y2DFLAXJ2DF2AXYIDFIBAJIBXSIGHACKXDEL TA2KX2 
+Y2DFIAXJ2DFCAXVY1IBXJIDFIBXSIGHA2XX2XSIGMAIX 
DELTA2*X2-Y2DFIAXJ2ZDFZAKYIBXJIDFIBXSIGMA2ZXX2*% 
DELTA2X¥2-Y2DFIAXJ2DF2ZAKY1BXJIDFIBXSIGMA2* 
SIGMALXDELTA2*®2+Y2DFIAXJ2DF2AXY1BXJIDFIB* 
SIGHAZ*XDELTAZ¥¥2+Y2DFIAXIJ2DEIE<yY LDR ZA 
SIGMACK*2*SIGMALXDELTA2XDELTAI-Y2DFIAXJ2DFIBXYIDF2A 
XJLBXSIGMAZX®2XDELTA2XDELTA1L+ANS8 

“ANS6=Y2DF24¥J2DF1AXY1 DF1BXJ1BXSIGMA2*X2XSIGMAL 
X®DELTA2kX2-Y2DF2ZAXJ2DFIAXYIDFIBXJLBXSIGMAZRR2 
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MOQ ANAND AND ANNA ANAND 


XDELTAZXX2-Y2DFZAXJ2DFLAXY1DFIBXJ1BXSIGMA2x 
SIGMALXDELTA2XX2+Y2DFZAXJ2DFIAXYIDFIBXJIBX 
SIGMAZXDELTA2XX2-Ye2DFZAXJZDFIAXYIBXJIDFIBX 
SIGMAZKAZXSIGMALXDELTAZ*X2+Y2DF2ZAXJZDFIAXY1 BX 
JIDFIBXSIGMACKX2XDELTA2ZXX2+Y2DFZAXJZDFLAXY1BxX 
JIDFIBXSIGMA2XSIGMALXDELTA2*X2-Y2DFZAXJ2ZDFIAXY1B 
XJIDFIBXSIGHAZXDEL TAZXX2-Y2DFZAXJ2ZDFIBXYI1DFIA 
XJLBXSIGMAZCARZXSIGMALXDE: TAZXX24+2.xXY2DFZAXJ2DF1B 
XYLDFLAXJIBXSIGMAZXSIGMA1LXDELTA2*X2-Y2DF2Ax 
JZ2ZDFIBXYIDFLAXJIBXSIGMALADELTAZXX2+Y2DF2AxX 
JZ2ZDFIBXYIBXJLDFIAXSIGMAZHX2XSIGMAILXDELTA2*XX2-2 .x 


_ Y2DFZAXJ2DF1BXY1BXJ1DF1LAXSIGMA2*SIGMA1XDELTA2 


XXZ+YZDFZAXJZDFIBXYLBXJIDFIAXSIGMA1XDEL TAZKX2 
+Y¥Y2DFZAXJ2BRYIDFLAXJIDFIBXSIGMA2XSIGMA1 XX2x 
DELTA2XDELTAIL-Y2DF2AXJ2ZBXY1DFLAXJ1DFIBXSIGMAZ2 
XSIGMALXDELTA2XDELTAL+ANS7 


“ANS5=AXOMEGAXANS6X(0.0,-1.0) 


IMPEDANCE 

ZZ = ANSI/ANSS5 

REAL PART OF IMPEDANCE 

ZR = REALCZZ) 

IMAGINARY PART OF IMPEDANCE 
ZI = AIMAGCZZ) 

MAGNITUDE OF IMPEDANCE 

Z = CABSCZZ) 

DB OF IMPEDANCE 

2UBe= 20.0xL0GIOCZ) 


RETURN 
END 
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APPENDIX C. 
OTHER FIGURES 
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Figure 66. Experimental Real Part of the Driving Point Impedance of 
a 6 in Radius Elastic Plate, With and Without Shaker and 
Mount Correction, in the Frequency Range of 5 to 105 Hz. 
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Figure 67. Experimental Imaginary Part of the Driving Point 
Impedance of a 6 in Radius Elastic Plate, With and 
Without Shaker and Mount Correction, in the Frequency 


Range of 5 to 105 Hz. 
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Figure 68. Experimental Real Part of the Driving Point Impedance of 
a 6 in Radius Viscoelastic Plate, With and Without Shaker 
and Mount Correction, in the Frequency Range of 100 to 
2000 Hz at a Temperature of 100.5 Deg. F. 
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Figure 69. Experimental Imaginary Part of the Driving Point 
Impedance of a 6 in Radius Viscoelastic Plate, With and 
Without Shaker and Mount Correction, in the Frequency 
sages of 100 to 2000 Hz at a Temperature of 100.5 
eg. F. 
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Figure 70. Effect of Shear Coefficient on the Theoretical Real Part of 
the Driving Point Impedance for a 5 in Radius Elastic 
Plate in the Frequency Range of 10 to 105 Hz. 
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Figure 71. Effect of Shear Coefficient on the Theoretical Imaginary 
Part of the Driving Point Impedance for a 5 in Radius 


Elastic Plate in the Frequency Range of 10 to 105 Hz. 
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Figure 72. Effect of Shear Coefficient on the Theoretical Real Part of 
the Driving Point Impedance for a 5 in Radius Elastic 
Plate in the Frequency Range of 100 to 2000 Hz. 
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Figure 73. Effect of Shear Coefficient on the Theoretical Imaginary 
Part of the Driving Point Impedance for a5 in Radius 
Elastic Plate in the Frequency Range of 100 to 2000 Hz. 
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‘igure 74. Effect of Shear Coefficient and Material Characteristics on 
the Theoretical Real Part of the Driving Point Impedance 
for a 5 in Radius Viscoelastic Plate in the Frequency 
Range of 10 to 105 Hz at a Temperature of 80.0 Deg. F. 
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Figure 75. Effect of Shear Coefficient and Material Characteristics on 
the Theoretical Imaginary Part of the Driving Point 
Impedance for a 5 in Radius Viscoelastic Plate in the 
Frequency Range of 10 to 105 Hz at a Temperature of 80.0 
Deg. F. 
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Figure 76. Effect of Shear Coefficient and Material Characteristics on 
the Theoretical Real Part of the Driving Point Impedance 
for a 5 in Radius Viscoelastic Plate in the Frequency 
Range of 100 to 2000 Hz at a Temperature of 77.3 Deg. F. 





140 


2000.0 


Pe HeH Oe eee ewer se ese eeesrMeseeseeeseet Min gperueeMeesesasereee 


e 
Com emer eee renee sere ea eetrseeteeeese ses eeseseasesemesesshbeeretes 


1 
FREQUENCY (IIZ) 





© 

S 

ee | ee co ot PU EL | che = pce cena re mrctctere ia ce sd a8 voica wien’ ciavainiovevaraya eis e.siela © ovals etsloeie ore cisdisiaieiSeieieicie icicle oS 

@ 

! o 

| S 

MI NN ae eon ene vo veces teeta gasses ne see+s sae 227) ad 

: © 
= : 
— ' ; 
ae 
- -— ew 
—_ “A+ 1rF 1 n 

= aw: pod | eee - = 

S adn s = 

oa tm ' (ped » = 

(me <r a = 
‘or I. 
poe | a 
ome | is 
caeoeayen 

nw =o | 

ome * rr So 

O ee ee eee = 
be ba Reelin | 

Bae a. N 
— |, — a 
= Cee ae 
— 

o 

ca) 

S 

= 

j . j i 
0°Se 0O°cl O-o- O°cT- 0o'-ce= 


0's 
(NI/OUS-AGT) JONVAIdNI 
Figure 77. Effect of Shear Coefficient and Material Characteristics on 
the Theoretical Imaginary Part of the Driving Point 
Impedance for a 5 in Radius Viscoelastic Plate in the 
Frequency Range of 100 to 2000 Hz at a Temperature of 
77.3 Deg. F. 
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Figure 78. Experimental Real Part of the Driving Point Impedance of 
a 5 in Radius Elastic Plate, With and Without Shaker and 


Mount Correction, in the Frequency Range of 5 to 105 Hz. 
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Figure 79. Experimental Imaginary Part of the Driving Point 
Impedance of a 5 in Radius Elastic Plate, With and 
Without Shaker and Mount Correction, in the Frequency 
Range of 5 to 105 Hz. 
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Figure 80. Experimental Real Part of the Driving Point Impedance of 
a 5in Radius Viscoelastic Plate, With and Without Shaker 
and Mount Correction, in the Frequency Range of 5 to 105 
Hz at a Temperature of 80.0 Deg. F. 
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Figure 81. Experimental Imaginary Part of the Driving Point 
Impedance of a 5 in Radius Viscoelastic Plate, With and 
Without Shaker and Mount Correction, in the Frequency 
Range of 5 to 105 Hz at a Temperature of 80.0 Deg. F. 
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Figure 82. Experimental Real Part of the Driving Point Impedance of 
a 5in Radius Viscoelastic Plate, With and Without Shaker 
and Mount Correction, in the Frequency Range of 100 to 
2000 Hz at a Temperature of 77.3 Deg. F. 
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Figure 83. Experimental Imaginary Part of the Driving Point 
Impedance of a 5 in Radius Viscoelastic Plate, With and 
Without Shaker and Mount Correction, in the Frequency 
Range of 100 to 2000 Hz at a Temperature of 77.3 Deg. F. 
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Figure 85. Theoretical and Experimental Comparison of the 
Imaginary Part of the Driving Point Impedance for a 5 in 
Radius Elastic Plate in the Frequency Range of 5 to 
105 Hz. 
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Figure 86. Theoretical and Experimental Comparison of the Real Part 
of the Driving Point Impedance for a 5 in Radius Elastic 
Plate in the Frequency Range of 100 to 2000 Hz. 
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Figure 87. Theoretical and Experimental Comparison of the 
Imaginary Part of the Driving Point Impedance for a 5 in 
Radius Elastic Plate in the Frequency Range of 100 to 
2000 Hz. 
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Figure 88. Theoretical and Experimental Comparison of the Real Part 
of the Driving Point Impedance for a 5 in Radius 
Viscoelastic Piate in the Frequency Range of 5 to 105 Hz at 
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Figure 89. Theoretical and Experimental Comparison of the 
Imaginary Part of the Driving Point Impedance for a 5 in 
Radius Viscoelastic Plate in the Frequency Range of 5 to 
105 Hz at a Temperature of 80.0 Deg. F. 
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Figure 90. Theoretical and Experimental Comparison of the Real Part 
of the Driving Point Impedance for a 5 in Radius 
Viscoelastic Plate in the Frequency Range of 100 to 2000 Hz 
at a Temperature of 77.3 Deg. F. 
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Figure 92. Effect of Shear Coefficient on the Theoretical Real Part of 
the Driving Point Impedance for a 4 in Radius Elastic 
Plate in the Frequency Range of 10 to 105 Hz. 
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Figure 93. Effect of Shear Coefficient on the Theoretical Imaginary 
Part of the Driving Point Impedance for a 4 in Radius 
Elastic Plate in the Frequency Range of 10 to 105 Hz. 
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Figure 95. Effect of Shear Coefficient on the Theoretical Imaginary 
Part of the Driving Point Impedance for a 4 in Radius 
Elastic Plate in the Frequency Range of 100 to 2000 Hz. 
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Figure 96. Effect of Shear Coefficient and Material Characteristics on 
the Theoretical Real Part of the Driving Point Impedance 
for a 4in Radius Viscoelastic Plate in the Frequency 
Range of 10 to 105 Hz at a Temperature of 78.5 Deg. F. 
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Figure 97. Effect of Shear Coefficient and Material Characteristics on 
the Theoretical Imaginary Part of the Driving Point 
Impedance for a 4 in Radius Viscoelastic Plate in the 
Frequency Range of 10 to 105 Hz at a Temperature of 
78.5 Deg. F. 
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Figure 98. Effect of Shear Coefficient and Material Characteristics on 
the Theoretical Real Part of the Driving Point Impedance 
for a 4 in Radius Viscoelastic Plate in the Frequency 
Range of 100 to 2000 Hz at a Temperature of 74.5 Deg. F. 
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Figure 99. Effect of Shear Coefficient and Material Characteristics on 
the Theoretical Imaginary Part of the Driving Point 
Impedance for a 4 in Radius Viscoelastic Plate in the 
Frequency Range of 100 to 2000 Hz at a Temperature of 
74.5 Deg. F. 
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Figure 100. Experimental Real Part of the Driving Point Impedance of 
a 4 in Radius Elastic Plate, With and Without Shaker and 
Mount Correction, in the Frequency Range of 5 to 105 Hz. 
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Figure 10L. Experimental Imaginary Part of the Driving Point 
Impedance of a 4 in Radius Elastic Plate, With and 
Without Shaker and Mount Correction, in the Frequency 
Range of 5 to 105 Hz. 
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Figure 102. Experimental Real Part of the Driving Point Impedance of 
a 4 in Radius Viscoelastic Plate, With and Without Shaker 
and Mount Correction, in the Frequency Range of 5 to 105 
Hz at a Temperature of 78.5 Deg. F. 


166 


105.0 


FREQUENCY (IZ) 


9.0 





Oe 0°Z O71 0°0 
(NI/O3S—AG1) JONVAIdNI 


Figure 103. Experimental Imaginary Part of the Driving Point 
Impedance of a 4 in Radius Viscoelastic Plate, With and 
Without Shaker and Mount Correction, in the Frequency 
Range of 5 to 105 Hz at a Temperature of 78.5 Deg. F. 
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Figure 104. Experimental Real Part of the Driving Point Impedance of 
a 4in Radius Viscoelastic Plate, With and Without Shaker 
and Mount Correction, in the Frequency Range of 100 to 
2000 Hz at a Temperature of 74.5 Deg. F. 


168 


2000.0 


143 


1240.0 
CY (IZ) 


1050.0 


860.0 
FREQUI 


670.0 


480.0 


0.0 


ao 





100.0 





0°01 O's 0°0 o’c- 0°0T- 
(NI/OAS-AGT) JONVAAdNI 
Figure 105. Experimental Imaginary Part of the Driving Point 
Impedance of a4 in Radius Viscoelastic Plate, With and 
Without Shaker and Mount Correction, in the Frequency 
Range of 100 to 2000 Hz at a Temperature of 74.5 Deg. F. 
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Figure 106. Theoretical and Experimental Comp 
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Figure 107. Theoretical and Experimental Comparison of the 
Imaginary Part of the Driving Point Impedance for a 4 in 
Radius Elastic Plate in the Frequency Range of 5 to 
105 Hz. 
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Figure 108. Theoretical and Experimental Comparison of the Real Part 


of the Driving Point Impedance for a 4 in Radius Elastic 
Plate in the Frequency Range of 100 to 2000 Hz. 
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Figure 109. Theoretical and Experimental Comparison of the 
Imaginary Part of the Driving Point Impedance for a 4 in 
Radius Elastic Plate in the Frequency Range of 100 to 
2000 Hz. 
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Figure 110. Theoretical and Experimental Comparison of the Real Part 
of the Driving Point Impedance for a 4 in Radius 
Viscoelastic Plate in the Frequency Range of 5 to 105 Hz at 
a Temperature of 78.5 Deg. F. 
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Figure 111. Theoretical and Experimental Comparison of the 
Imaginary Part of the Driving Point Impedance for a 4 in 
Radius Viscoelastic Plate in the Frequency Range of 5 to 
105 Hz at a Temperature of 78.5 Deg. F. 
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